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FOREWORD 


This repoft describes the results of work performed by Mechanical Technology 
Incorporated (MTI), Latham, New York, under Task I of NASA Contract NAS3-15709. 
Experimental data reduced under this Task was obtained and reported in NASA 
Contractor Report CR-2177, October 1971, under Contract NASw-1713. 

The NASA Project Manager for this investigation was Mr. James Dunn of the Brayton 
Cycle Branch of the Lewis Research Center, Cleveland, Ohio. The MTI Project 
Manager was Dr. Robert H. Badgley. Mr. Juergen M. Tessarzik of MTI was responsi- 
ble for the experimental portion of the investigation, while Dr. Thomas Chiang 
was responsible for the analytical portion. 

The computer program used to calculate axial random vibration response probabil- 
ities for the BRU simulator rotor-bearing system was developed by Dr. Thomas 
Chiang . 

The authors are especially grateful to Mr. Walter Spodnewski of MTI for his 
skillful assistance in performing the power spectral density and probability 
laboratory analyses of the experimental test data. 
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ABSTRACT 


The random vibration response of a gas-bearing rotor-support system has been 
experimentally and analytically investigated in the amplitude and frequency 
domains. 

The NASA Brayton Rotating Unit (BR.U), a 36,000 rpm, 10 KWe turbogenerator had 
previously been subjected in the laboratory to external random vibrations, and 
the response data recorded on magnetic tape. This data has now been experimen- 
tally analyzed for amplitude distribution and frequency content. The results 
of the power spectral density analysis indicate strong vibration responses for 
the major rotor-bearing system components at frequencies which correspond 
closely to ,their resonant frequencies obtained under periodic vibration testing. 
The results of amplitude analysis indicate an increasing shift towards non- 
Gaussian distributions as the input level of external vibrations is raised. 

Analysis of axial random vibration response of the BRU was performed by using a 
linear three-mass model. Power spectral densities, the root-mean-square value 
of the thrust bearing gas-film thickness variation, and the probability of 
thrust bearing surface contact were calculated for specified input random 
excitation . 
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SUMMARY 


The axial random vibration response of the NASA Brayton Rotating Unit (BRU) was 
analyzed with the BRU represented by a three-mass linear system. A computer 
program was extended to permit calculation of power spectral densities, the root- 
mean-square value of the thrust bearing gas -film thickness variation and the 
probability of thrust bearing surface contact for specified input random excita- 
tion. 

Experimental analysis of the random BRU simulator response data has also been 
accomplished. Probability density and distribution functions have been computed 
for the data for all representative rotor-bearing system displacement sensors at 
low and high external vibration input levels. Similarly, a frequency analysis 
(displacement power spectral density) for the same data has been obtained. 

The probability analysis has revealed a wide range of BRU system component 
responses to the externally-applied random vibration environment. Some compo- 
nents reacted with normal (Gaussian) or nearly-normal responses at all vibration 
test input levels, while others exhibited drastic changes towards non-Gaussian 
responses for increasing vibration input levels. Generally, implications about 
the physical behavior of bearing system parts can be drawn from the examination 
of their amplitude response analysis. 

The frequency analysis of the BRU simulator response data to external random 
vibration inputs has revealed a general similarity in the response patterns 
obtained for random and sinusoidal vibration inputs. Major system resonances 
occurred in the frequency range between 60 and 400 Hz. The gain factor (ratio 
of vibration response level to input level) for most system components under 
random vibration was well under unity, with the rotor axial vibration generally 
exhibiting the highest relative response. 

Evaluation of the BRU simulator bearing system component responses to sinusoidal 
vibration inputs over the frequency range 5 to 2000 Hz indicated a large number 
of fairly closely spaced resonances in the frequency range 60 to 400 Hz. These 



appear to be due to resonant behavior of the rotor mass supported on the bearing 
gas films . 

The major factors contributing to the multiplicity of observed 'resonances appear 
to include unequal journal bearing gas-film stiffnesses in the direction through 
and at right angles to the pivot of the flexure-supported pad, and unequal stiff- 
nesses of the two journal bearings. Major lateral resonances appear to be of the 
rigid body type common to two-bearing rotor systems. Although the absolute dis- 
placement response level varies for individual parts (depending upon the direction 
of applied vibration relative to the direction of freedom of the particular part^ 
no single resonance was found to be especially severe. 



-3 


INTRODUCTION 

The adaption of Brayton Cycle Space Power Systems to environments involving 
externally- imposed vibrations must logically begin with an evaluation of the 
capabilities of existing equipment. 

The Brayton Rotating Unit (BRU), a gas turbine - driven alternator of approximately 
10 KWe maximum output with a proven steady-state performance, had previously been 
subjected to sinusoidal and random vibration environmental test conditions (Ref- 
erence 1). Test results appeared encouraging, but were not fully evaluated, 
especially in the area of random vibration. 

In this report, previously recorded data from vibration tests of the BRU are sum- 
marized and evaluated for the case of sinusoidal vibration input conditions. 

Data obtained for random vibrational input conditions has been experimentally 
analyzed extensively in the frequency and amplitude domains (displacement power 
spectral density and probability density, respectively).. 

A computer program previously developed (Reference 1) for the analytical predic- 
tion of the axial rotor response in terms of acceleration power spectral density, 
under external random vibrations, has been extended to permit the calculation of 
the probabilities that thrust bearing gas- film thickness variations will be 
greater than specified values (probability density) for a specified design and 
vibration input condition. 

Previously calculated acceleration power spectral density results are compared 
with the newly analyzed experimental data. Vibration excursion probabilities 
have been calculated for the rotor-thrust bearing combination under external 
random vibrations in the rotor-axial direction. 



PROBABILITY OF THRUST .BEARING. SURFACE CONTACT DUE 
TO:. EXTERNAL. AXIAL RANDOM VIBRATION. 


Brayton-cycle turbomachinery may be subjected to externally-imposed random vibra- 
tion under general usage conditions. Under random excitation, the internal 
machine elements, and particularly the gas bearings, will respond dynamically. 

An analysis of the axial dynamic response of the Brayton Rotating Unit (BRU) due 
to external random excitation was presented in Reference 1. For any prescribed 
random input, the response in terms of the power spectral density and the rms 
(root-mean-square) value of the thrust bearing gas-film thickness variation was 
obtained. In this report, the analysis is extended to include a probability 
analysis of the thrust bearing gas-film thickness variation. 

Theoretical Formulation of Probability Analysis 

Consider a random time-history record y(t) . The probability density function of 
y(t) is defined as the probability that y(t) will fall within a specific range at 
any given: instant of time. In a time duration T, let Tiy be the total amount of 
-.time that y(t) falls within the range y and y + Ay . Then, as T approaches in- 
finity, the ratio T x /T will approach the probability that y(t) falls inside the 
range (y, y + Ay). In mathematical terms, 

1 im J~T ~j 

Prob[y < y(t) < y + Ay] = ^ |jj?j (1) 


The probability density function is defined as 

n / v N _ l.im Prob[y < y(t) < y + Ay] 

PKy r Ay->0 Ay 

( 2 ) 

= Frob [y < y(t) < y + dy] 

dy 


This probability density function is to be distinguished from the probability 
distribution function P(y) which is defined in the following; 
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Thus , the probability of y being less than -» is zero and the probability of y 
being less than “ is unity. 



In the probability analysis of the BRU thrust bearing gas-film thickness varia- 
tion due to random excitation, the primary interest is centered on the probability 
that the film thickness will be zero, which is of course the probability of bear- 
ing surface contact, P' . Let "h" be the mean gas-film thickness. Then, from 
Equation (3) 

P(h) = Prob[y(t) < h] = ^ p(y , )dy' (5) 

-00 

Thus, the probability of bearing surface contact is, 

P' = Prob[y(t) > h] 

CO 

P(y')dy' 

= 1 ~ P(h) 

In Appendix A the analysis of the response power spectral density and rms value, 
which was shown in Appendix C of Reference 1, is recapitulated for ready reference. 




The analytical model used to study the BRU axial random vibration problem is a 
linear constant-parameter system. This important assumption permitted the re- 
sponse analysis to be carried out in a simple manner. The linearity asumption 
is also essential for the probability analysis. 

I 

The input random excitation applied by the vibration' system to the BRU during 
testing has a normal or Gaussian, distribution. The response of a linear system 

must, therefore, also have a Gaussian distribution. Based on this distribution, 

. J 

the probability of thrust bearing surface contact may be calculated. A detailed 
description of the probability analysis is presented in Appendix A. 

Appendix B contains a complete listing of the computer program which implements 
the analysis of Appendix A. The values of the masses, stiffnesses, and damping 
coefficients are inputs to the computer program. Also included in Appendix B 
are instructions for preparing input data cards, sample input, and sample 
printed output. 

Results of Probability Calculation of BRU Thrust Bearing Surface 
Contact Due to Random Excitation 

The foregoing response and probability analyses were conducted for the case of 
axial random vibration applied to the Brayton Rotating Unit with shaft rotation. 
Before giving the results of the analysis of the probability of thrust bearing 
surface contact, it is appropriate to describe the input to the program. 
(Nomenclature refers to Figure 77.) 

Inputs 

m^ = mass of the BRU housing structure = 96 lb. 

m 2 = mass of the gimbal and thrust bearing = 7 lb . 

m^ = mass of the rotor = 21 lb. 

k^ = stiffness of connector between the BRU and shake table = 300,000 

lb/in. 

k 2 = stiffness of the gimbal support = 250,000 lb/in. (from page 79 

of Reference 2) 

k^ = stiffness of the gas thrust bearing = 113,500 lb/in. 
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c = damping coefficient of connector = 0.02 lb sec/in. 

= damping coefficient of gimbal support = 0.02 lb sec/in. 

c^ = damping coefficient of gas thrust bearing = 84 lb sec/in. 


The above values of the masses, stiffnesses and damping coefficients are the 
nominal values and are summarized in Table I. 




TABLE I 


Nominal Values of the Lumped Masses, Stiffnesses, 
and Damping Coefficients of an Axial Vibration 
Model of the Brayton Rotating Unit 

Element 

m (lb) 

k (lb/in.) 

c (lb sec /in.) 

1 

96 

300,000 

0.02 

2 

7 

250,000 

0.02 

3 

21 

113,500 

84.00 


It should be noted that the values of k^, and c^ are only nominal values. 

With respect to the stiffness k^ , there are three connectors between the BRU 
and the shake table. Each is bolted down at one end to the shake table. When 
the unbolted end of the connector tends to move up, the connector behaves like 
a cantilever with a calculated stiffness of 75,000 lb/in. at the BRU mount 
point. But for downward motion the cantilever arm is considerably shorter and 
thus much stiffer. This stiffness is calculated to be about 1,200,000 lb/in. 

For the purpose of this analysis, a nominal value of k^ = 300,000 lb/in. is 
chosen. 

The stiffness and damping characteristics of the gas thrust bearing are functions 
of frequency and bearing clearance. They have been calculated by an existing 
gas-bearing computer program and tabulated in Tables II and III. It is seen 
that they both vary widely with frequency and bearing clearance. Since the 
random vibration analysis is limited to a linear, time-invariant system, the 
stiffness and damping coefficients must have fixed values. The nominal values 
of k 3 and c^ are chosen at a frequency of 100 Hz and with a film thickness of 





TABLE II 

Dynamic Stiffness (lb/in.) of Double-Acting 
Gas Thrust Bearing 


NMin.) 

o>(Hz) 

.0002 

.0005 

.00105 

100 

477,000 

86,900 

18,900 

500 

523,000 

93,700 

20,000 

1,000 

652,000 

113,500 

22,700 

2,000 

950,000 

166,000 

33,600 


TABLE III 

Dynamic Damping Coefficient (lb sec /in.) 
of Double-Acting Gas Thrust Bearing 


^S^hCin.) 

o>(Hz)^'' v 'v v ^ 

.0002 

.0005 

.00105 

100 

413 

92 

33 

500 

392 i 

90 

33 

1,000 

338 

84 

33 

2,000 

191 

66 

32 
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0.0005 inch (the total axial clearance in the thrust bearing is 0.0021 inch). 

The test objective power spectral density S Q is shown in Figure 1. During actual 
testing of the BRU, to be described below, it was found that frequent touching 
of the gas thrust bearing surfaces (possibly to the point of failure) would have 
occurred if the attempt had been made to subject the BRU to the full level of 
this vibration specification. Tests were therefore performed with reduced input 
PSD values (i.e., the shape of the PSD was preserved, but the magnitude of the 
PSD was reduced by a constant factor, a, across the whole spectrum). From the 
definitions of the power spectral density and the root-mean-square value, the 
root-mean-square value is reduced by a factor of Ja when the power spectral 
density is reduced by a factor of a across the whole spectrum. 

The root-mean-square (rms) acceleration of the power spectral density in Figure 
1 is 5.4g. Testing was performed on the BRU at reduced PSD with the same shape, 
but with several rms acceleration levels including 0.54g. The response PSD is 
calculated based on this 0.54g rms acceleration. Thus, ,Ja = ~ ^ = 0.1 and 
O' = 0.01. 

With the above input values of masses, stiffnesses and damping coefficients, and 
input PSD of Figure 1 with a = 0.01, the response PSD’s Sy^, Sy 2 and Sy^ were 
calculated. Note that y^, y 2 and y^ are the relative displacements between m^ 
and base, m 2 and m^, and m 2 and m 2 , respectively. The area underneath the PSD 
curve represents the mean square value of the respective quantity. 

2 -8 2 

E(y^) = mean square value of y^ = 9.125 x 10 in. 

2 -82 

E"^ 2 ^ = mean square value of y^ = 1.753 x 10 in. 

2 -8 2 

E(y^) = mean square value of y^ = 5.164 x 10 in. 

The square root of the mean square value is the root-mean-square (rms) value. 

Thus , 

-3 

rms of y^ = 0.226 x 10 in. 

-3 

This is to be compared with the measured value of 0.35 x 10 inch. 

The rms of y 2 is the statistical average of the thrust bearing gas-film thickness 



variation. On the other hand, the mean value of is zero if the response to 

random excitation has already reached a steady state. Therefore, the rms of y^ 

is also the standard deviation of which is denoted by a,, . 

'3 ^3 


If the input random excitation is low, the system remains linear. As stated 
earlier, the random input from the shake table to the BRU has a normal or Gaussian 
distribution, and for a linear system, the response to an input with Gaussian 
distribution also has a Gaussian distribution. Thus, knowing the rms value of 
y^ , (which in this case is also the standard deviation of y^) and that y^ has a 
Gaussian distribution, the probability of thrust bearing surface contact may be 
calculated and tabulated as a function of mean gas-film thickness. 


Mean Gas-Film Probability of Thrust Bearing 

Thickness (in.) Surface Contact 


0.0005 

0.0006 

0.0007 


0.0139 

0.0042 

0.0011 
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INSTRUMEN TA-TTON for expe rimental -analysis 
T': OF ..RANDOM VIBRATION! TEST . MIA’ . 

Introduction 

The BRU simulator had previously been subjected to external sinusoidal and 
random vibrations which were applied in the rotor axial direction and in a 
direction perpendicular to the rotor axis . The response of the rotor and of 
many bearing system components was recorded on magnetic tape for later analysis 
(Reference 1) . The experimental analysis presented in this report draws upon 
the data obtained and cataloged in Reference 1. A complete description of the 
Test Apparatus (BRU simulator) , the test instrumentation and a description of 
vibration tests performed may also be found in Reference 1„ For convenience, 
an abbreviated description of the test apparatus and the instrumentation is 
repeated here. 

Test Apparatus 

A schematic of the BRU simulator with the rotor in the vertical orientation is 
shown in Figure 2. The rotor, from the top down, consists of the following 
components ; 

1. The (simulated) compressor wheel, containing the cold gas simulator 
drive turbine . 

2. The thrust runner. 

3. The compressor end journal. 

4. A center section (homopolar generator rotor). 

5. The turbine end journal. 

6. The (simulated) gas turbine wheel. 

The three-pad configuration of the gas film journal bearings is shown in cross- 
section in Figure 2. The pads are individually supported by pivots, with one of 
the pivots in each bearing assembly mounted on a mechanical flexure. The flexure, 
which has a radial stiffness of approximately 2000 pounds per inch, permits the 
bearing to accommodate radial centrifugal growth of the journal and differential, 
radial, thermal expansion between the various bearing parts. 



The double-acting thrust bearing is supported by a gimbal assembly with two 
sets of gimbals 90° apart. These permit the thrust bearing to align itself with 
the thrust runner of the rotor. The surface geometry of the thrust plates con- 
sists of nine equal sectors of 39° each with narrow radial grooves separating the 
sectors. Each sector in turn consists of a slightly depressed sector of 15° arc 
followed by the raised part in the direction of rotor rotation. Design details 
for the complete journal and thrust bearing assemblies are given in Reference 2. 

The journal bearings, as well as both thrust bearing plates, have hydrodynamic 
and hydrostatic operating capabilities. For hydrostatic operation, each journal' 
bearing pad and each raised thrust plate sector is equipped with a small gas 
supply hole, connected to an outside gas source. Hydrostatic bearing operation 
is normally employed to separate the bearing and rotor surfaces at start-up 
until rotor speed is high enough to produce a sufficiently large gas film for 
safe bearing operation. At the design speed of 36,000 rpm and above, the rotor 
operates hydrodynamically . The BRU design also requires hydrostatic bearing 
operation at shutdown to prevent rubbing contact between the rotor and the 
bearing surfaces at speeds too low to provide hydrodynamic film separation of 
the bearing surfaces . Prior to the vibration tests described in this report, 
all bearing and corresponding rotor surfaces had been coated by MTI with chrome 
oxide and refinished to the original dimensional specifications. This was done 
to improve the tolerance of the bearing surfaces against accidental rubs which 
might occur in the course of the vibration tests . 

Figure 3 is a photograph of the BRU simulator in the vertical support fixture 
mounted on the vibration table ready for testing. The simulator is held in the 
support fixture without isolators by three mounting brackets extending radially 
outward from the simulator casing and resting on the rim of the vertical support 
fixture. The flexible air hoses visable in Figure 3 serve to pressurize sepa- 
rately each of the journal bearings, the thrust bearing plates, the bearing 
cavity, the thrust loader chamber and the cold gas drive turbine. Nitrogen at 
room temperature was used for hydrostatic bearing operation and to provide the 
ambient gas atmosphere during all hydrodynamic tests. 

Figure 4 .is a photograph of the BRU simulator in the transverse test fixture 
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which is attached to the shake table surface. The shake table has been rotated 
into a horizontal position so that the direction of vibration is perpendicular 
to the vertically oriented axis of rotation of the BRU. The same three mounting 
brackets without isolators are used for connecting the BRU simulator to the 
transverse test fixture as were used on the vertical fixture. 

A small-scale schematic of the BRU simulator in the transverse test fixture is 
shown on all data graphs obtained for this arrangement. As the schematic in- 
dicates, the direction of applied vibration is in line with the instrumentation 
turret and thus passes also through the pivot of the flexibly-mounted pad in 
each journal bearing. 

Test Instrumentation 

The dynamic response measurements on the rotor-bearing system were made with 
capacitance probes. These probes may be divided for convenience into two cat- 
egories; those measuring displacements of the rotor or some bearing component 
relative to the simulator casing, and those measuring gas-film thickness varia- 
tions between the rotor and the thrust bearing or between the rotor and any of 
the journal bearing pads. From a total of four probes measuring thrust bearing 
film thickness variations, only two data signals from two probes 90° apart cir- 
cumferentially have generally been recorded. These two probes were axially in 
line with probes recording motions of the gimbaled thrust bearing relative to 
the casing. 

Gas-film thickness variations between the rotor and individual journal bearing 
pads, were measured with capacitance probes built into each pad. These probes, 
which were installed by MTI prior to tests, were located next to the pivot on 
the inboard side of each bearing pad (a total of six). 

Two sets of two probes were used to measure rotor displacement relative to the 
simulator casing at two locations inboard of the journal bearings. (The angular 
orientation of these capacitance probes is indicated in Figure 5.) The radial 
motions of the two. flexures supporting each one of the three pads of each journal 
bearing were measured near the pivot location relative to the casing. Pad 
motions at the pad leading edge were recorded for one solid-mounted and one 



flexure-mounted pad, A schematic of the probe locations in relation to the 
rotor is shown in Figure 5. 

Temperatures from nine locations, including both flexure-mounted journal pads, 
one solid-mounted pad at the compressor-end journal, and one thrust plate were 
recorded by automatic print-out throughout the tests. 

Analyzers 

The random vibration response amplitudes, which had been recorded on magnetic 
tape during earlier tests, were analyzed in three steps, each employing a 
distinctly different type of analyzer. 

The root-mean-square values of the random displacement signals were extracted 
through application of a true rms reading voltmeter. The model used was a 
Hewlett-Packard 7562A log voltmeter-converter. (Many of the simpler commercial 
voltmeters take advantage of the fixed relationship between peak or mean abso- 
lute values to the instantaneous values of a sine signal — which normally 
represents an a-c voltage. These instruments are much less expensive than true 
rms detection circuits, but will produce inaccurate readings if they are used to 
measure random signals.) 

The probability analyzer used on the amplitude analysis of the tape recorded 
BRU simulator response data was an instrument of the latest development . (Sig- 
nal Analysis Industries Corporation, Model SAICOR 43-A.) The all-digital high- 
speed processing instrument provided on-line, real-time computation which did 
not require the use of tape loops for data processing. The probability functions 
are determined at 400 points with the horizontal axis of the display representing 
voltage and the range and sensitivity being determined by the setting of the in- 
put attenuators. The averaging process can be accomplished linearly or in the 
digital exponential mode at a minimum sampling rate of At = 0.2 microsecond. 

The computed functions are displayed on an oscilloscope screen and may be plotted 
from instrument memory after completion of the computation. A typical photo of 
the oscilloscope display and the corresponding plot produced on an x-y plotter 
are shown in Figure 35 . - 
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The tape recorded BRU random data was further analyzed with the aid of a real- 
time frequency analyzer and its companion instrument, an ensemble averager. 
(Spectral Dynamics Corporation, SD 301 Real-Time Analyzer and SD 302 Ensemble 
Averager.) This powerful and highly efficient combination permitted the averaging 
of up to 256 real-time spectra of BRU data for frequency analysis with bandwidth 
resolutions between 0.3 and 6 Hz over frequency ranges between 100 and 2000 Hz. 
(These values are not the performance limits of the instrument, they rather 
represent the range of parameters used in the BRU data analysis and were largely 
dictated by the length of the available tape recorded random data. The applica- 
tion of bandwidth variation to the BRU data analysis is further explained and 
demonstrated at a later section of this report, titled "Frequency Analysis — 
Displacement Power Spectral Density.") 

The analysis results obtained by the Real-Time Analyzer are plotted on paper 
from instrument memory. 



RESULTS OF .EXPERIMENTAL ANALYSIS OF 
RANDOM VIBRATION TEST DATA 

Amplitude Versus Time Records 

Under random vibration conditions, visual inspection of the 'raw' displacement 
response data will yield only a limited amount of quantitative information. The 
amplitude limits may thus be determined, provided a long enough sample of data 
is being examined. A long data record is required to ensure that no gradual 
shift in the general signal amplitude level is taking place, even though the 
vibration input level is held constant. If such a change were detected, the 
data would be termed nonstationary and the experimental analysis methods used to 
obtain statistical data as presented here would not be applicable. 

Visual inspection of the raw data generated for the BRU simulator reveals that, 
within the maximum amplitude limits of the data the ratio of very large excur- 
sions to what appears to be the average signal level, varies from one data 
signal to the next. In fact, this ratio appears to vary drastically for the 
same signal when the level of the externally-applied vibration is changed. Some 
of these signals appear to be distributed so unevenly in level that it is quite 
safe to guess that their probability density function will be skewed — in other 
words, the occurrence of small and large amplitudes will not follow the bell- 
shaped "normal" or Gaussian distribution curve. 

An example of data of this type is shown in Figure 6. Both photographs represent 
the displacement response (relative to machine casing) of one of the f Ifexibly- 
mounted pivoted journal bearing pads. For the upper photograph the random vibra- 
tion input level was 0.54g (rms); and for the lower, the input level was 1.54g 
(rms) . For the low level input, the response was very nearly Gaussian, as in- 
dicated by the amplitude probability density curve shown next to it. (The 
probability curves shown in this figure are scale drawings obtained by experi- 
mental analysis from random data of which short segments are shown in the 
oscilloscope photographs.) When the vibration input level was raised (lower 
photo) the amplitude distribution shifted, with most amplitudes then occurring 
at a higher value and a few, very substantial excursions, at the low-value region 
(Physically, the low value excursions could either mean that the leading edge of 
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the pivoting pad came very close to hitting the rotating shaft, or conversely, 
that the shaft had moved far away — in the direction of the solid-mounted pad — 
and the :f lex-mounted pad was making a severe follow-up motion. In either case, 
the information from the photographs could be very valuable if vibration-induced 
pad instabilities were detected through shaft failure.) 

In some instances, only inspection of the time history of the displacement data 
will reveal otherwise unavailable system response information. A good example 
of this appears in the upper photograph in Figure 7, which shows the film thick- 
ness variations between the rotor and each of the flexibly-mounted bearing pads 
in each of the two journal bearings. The photo indicates that both pads behave 
in nearly identical fashion, executing motions that are comparatively smooth and 
do not change randomly in amplitude. It may also be noted that the amplitudes 
are out of phase, indicating apparently a "conical" vibration of the rotor. (To 
obtain the latter information, it was almost imperative to record both traces 
simultaneously. Due to the lack of periodicity in the signals, they must be 
recorded with phase and time integrity.) The second photo in Figure 7 records 
the much more typical response to random vibration of the film thicknesses be- 
tween the rotor and one of each of the two solid-mounted pivoted pads in each 
of the two journal bearings. 

In studies of rotating shafts, it is common to measure and display rotor motion 
in a plane perpendicular to the rotor centerline. While such motions are usually 
steady (not necessarily small) in a nonvibrating environment, they will be very 
unsteady (but not necessarily unstable) under random vibration conditions. In 
the left column in Figure 8 the unsteady or random nature of rotor motions is 
demonstrated through a series of photographs taken with identical film exposure 
times, which show the rotor moving through vastly different distances in identi- 
cal time spans. The right column in Figure 8 demonstrates, through photographs 
with increasingly longer exposure times, that it may be possible to observe the 
rotor vibrating relatively little during short and judiciously selected time 
periods. Longer and continuous observations will, however, show some of the 
very large excursions which, on a statistical basis, must be expected if they 
are not physically prevented in the system under test. 



An excellent example of distorted vibration response of a machine member due to 
limitations imposed upon its physical movement may be seen in Figure 9. Both 
photos show the gas film thickness variation between the thrust runner and the 
upper plate of the thrust bearing stator. Total thrust bearing clearance and 
the equilibrium operating position of the rotor is indicated to the left of each 
photo. (Differential gas pressures in the BRU simulator housing cause the thrust 
runner to assume an operating position which is above the axial centerline of 
the bearing, see Figure 10.) At low vibration input levels, the thrust runner 
vibrates essentially against the rather stiff gas film spring of the upper thrust 
plate without actually making meta 1- to-meta 1 contact. At the increased vibration 
level, the thrust runner moves the full distance between stator plates making 
frequent contact at both ends of its travel. (While thrust runner-to-stator 
contacts can only be inferred in Figure 9 from the fact that the thrust runner 
moves through a distance which equals as closely as can be determined the full 
thrust bearing clearance, temperature measurements in Figure 11 clearly prove 
that actual contact and intermittent rubbing was taking place.) 

Amplitude Analysis; Root-Mean-Square (rms) Values 

Calculation of the square root of the mean square (rms) values is the simplest 
form of data reduction on random signals. The procedure is well known (the rms 
value is equal to the standard deviation for quantities with a mean value of 
zero) and the value is readily determined for a random signal with a "true" rms 
reading voltmeter. Two sample plots (Figures 12 and 13) of random vibration 
amplitude rms values for vibration input values of 0.54, 0.75 and 1.52 g (rms) 
are reproduced here from Figures 83 and 86 of Reference 1. Similar plots were 
obtained for all measuring probes that were mounted inside the BRU simulator. 

The principal information that may be derived from the rms amplitude plots is 
the relationship between vibration input level (g-rms) and the amplitude re- 
sponse level (inch-rms) . Table IVa provides a summary of amplitude rms values 

at two vibration input levels [(0.54g (rms) and 1.52g (rms)] for all measurement 

probes with external vibrations applied to the BRU in the axial direction. The 
last column in Table IVa shows the ratio of rms amplitudes for the minimum (0.54g 
(rms)) and maximum (1.52g (rms)) vibration input levels. The amplitude response 
change is remarkably uniform, with the majority of the response ratios within 
20 percent of the input level ratio of 2.81. There are two exceptions for which 
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TABLE IV 

Root-Mean-Square Amplitudes at Two Random Vibration Input Levels 


(a) Axial Direction 


Probe i 
Number 



0 

0.52g (rms) 

1.52g (rms) 


4 

A 

u 

B 

„ 4 

B 

in. x 10 

in. x 10 

in o x 10 

A 

0.2 

0.375 

0.75 

2.0 

0.25 

0.375 

0.85 

2 .2 

0.3 

0.75 

2 .5 

3.3 

0.3 

1.25 

3.0 

2.4 

0.6 

1.2 

3.0 

2.4 

6.43 

1.1 

3.1 

2 .8 

0.68 

0.95 

2 .3 

2.4 

6.61 

1.1 

2.6 

2.4 

0.42 

1.1 

3.2 

2.9 

0.47 

3.6 

7.0 

1.9 

1.0 

2.0 

6.3 

3.1 

0.4 

3.25 

6.5 

2.0 

0.77 

2.0 

6 .6 

3.3 



(b) Transverse Direction 


Random Vibration Input Level 


0 

0,54g (rms) 

1.35g (rms) 


1.52g (rms) 

■ | 

Probe 

L 

A 

lx 

B 

u 

B 

C 

in. x 10 4 

s 

Number 

in. x 10 

in. x 10 

in. x 1.0 

A 

A 

A 

0.2 

1.1 

2 .5 

2.3 

2.8 

2 .5 

D 

0.29 

0.42 

0.8 

1.9 

N.A. 

N.A . 

C 

0.3 

0.75 

2 .2 

2.9 

2.75 

3.8 

F 

0.34 

2.0 

4.8 

2.4 

N.A. 

N.A. 

22 

0.72 

3.1 

9.0 

2.9 

11.0 

3.5 

21 

0.52 

2.3 

7.5 

3.3 

8.5 

. 3.7 

13 

0.8 

0.85 

1.8 

2.1 

2.1 

2 .5 

2 

0.75 

3.0 

8.0 

2.7 

9.0 

3.0 

4 

1.2 

3.5 

12 .0 

3.4 

N.A. 

N.A. 

20 

0.21 

1.9 

5.5 

2.9 

5.5 

2.9 

24 

0.61 

1.65 

3.45 

2.1 

N.A. 

N.A . 

17 

0.31 

1.5 

4.75 

2.3 

N.A. 

N.A. 

23 

1.05 

1.7 

3.75 

2.2 

N.A. 

N.A. 


































the increase in the vibration response is less than that of the input change: 
the thrust bearing probes (No. 17 and 20) and the probes in the two flexibly- 
mounted pivoted pads (A and D) . Both exceptions can be explained from known 
physical characteristics of the BRU, and from the respective amplitude time rec- 
ords (oscilloscope photographs, Figures 7 and 9). The flexibly-mounted pads 
have the capability to "track" the rotor to some extent, thereby sharing some of 
the rotor- induced vibrations with their respective support flexures. The dis- 
torted thrust bearing film thickness variations at higher vibration input levels 
were discussed earlier. Additional evidence for the view that at increased ran- 
dom vibration input levels the thrust runner will essentially bounce back and 
forth between the rotor plates (vibrate across the full thrust bearing clearance) 
is given in Table IV(b). Table IV(b) contains data similar to that in Table 
IV (a) except that it was obtained for the transverse direction of input vibra- 
tions. Since vibration response readings at the 1.54g (rms) input level were 
available for only a limited number of probes, additional values for the 1.35g 
(rms) vibration input have been included in this table. From this additional 
data, it is apparent that thrust runner rms amplitudes (probe 20) did not in- 
crease above the value reached at the input level of 1.35g (rms). The rms 
amplitude graph (Figure 139 of Reference 1) which includes test results for ran- 
dom vibration input levels as high as 1.70g (rms) confirms this observation. 

Table IV(b) also confirms the expectation that vibration inputs in the lateral 
direction will ellicit more severe gas film thickness variations in the solid- 
mounted journal bearings and larger vibrations of the flexure supporting the 
flexibly-supported pads than axial vibration inputs. 

Root-mean-square amplitude plots, as shown in Figures 12 and 13 have a very 
"noisy" appearance, because generally the damping built into the circuit of the 
measurement instrument is not large enough to eliminate the instantaneous fluc- 
tuations in the rms value. This could be considered as an advantage because the 
noisy looking rms curve also permits some quick, though superficial, assessment 
of changes in the nature of the random data signal. If the rms amplitude swings 
are small, the random vibration amplitude record could be expected to be rather 
uniform, which comparison of the records for 1.54g (rms) vibration input in 
Figures 13 and 9 will confirm. Small rms variations are, of course, not at all 



synonymous with small random vibration levels, as Figures 12 and 13 will show. 
Neither are the fluctuations of the rms amplitude curve an indication of the 
ratio between the random signal maximum amplitude envelope and the rms value. 

Theoretically, the random vibration maximum envelope should extend from — 00 to 
+ ® (Gaussian distribution), but shake table limitations (excursions are limited 
to + 3a about the mean) and physical restraints in the tested mechanism (BRU) do 
not allow a response that significantly exceeds + 3a. Once these restraints and 
the statistical uncertainty from observation of relatively short time records 
have been accepted, it is reasonable to question what the approximate ratio is 
between rms amplitude values and the random signal maximum envelope (amplitude- 
time record) for the set of response measurements at hand. Table V gives a 
listing of such values for the test case where external random vibrations were 
applied to BRU simulator in the axial direction. The average ratio of rms am- 
plitude value to maximum envelope in Table V is found to be five, or + 2.5 a. 
There are two very low values (of approximately + 1.5a) which were obtained for 
the two thrust bearing film thickness probes at th'e4higher (1.52g (rms)) vibra- 
tion input level. Of the four highest values, three were recorded for the hydro 
statically pressurized compressor-end journal bearing at the higher vibration 
input level and one, the highest of all (+ bo), was obtained for the turbine-end 
flexible pad-to-shaft film thickness variation. 

Amplitude Probability Analysis 

Determination of the probability density and distribution functions for a ran- 
dom signal is an important and necessary step in the data analysis progression. 
Knowledge of the type of amplitude distribution (normal or not normal) which 
flows from the analysis of the random data permits certain conclusions to be 
drawn relative to the physical phenomena underlying the recorded data. Moreover 
this knowledge can be used to determine some of the basic characteristics of the 
random data, which in turn may strongly influence the selection of procedures to 
be used for further analyses and interpretation of test results. 

The three most important basic characteristics of random data are the station- 
arity, the presence of periodicities and the normality. Stationarity is of con- 
cern because the analysis procedures used for stationary data are much simpler 
than those for nonstationary data, and do not apply at all to the latter. 
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Testing of random data stationarity should precede the amplitude probability 
analysis . For the case of the BRU simulator random data stationarity is conve- 
niently determined from the rms amplitude value plots discussed and referenced 
in the previous section. In general, the degree of stationarity as indicated by 
those plots, is quite satisfactory. For the vast majority of all recorded tests 
the rms response level remained constant at each external vibration input setting. 
Sample records were judged of sufficient length to permit differentiation of any 
nonstationary trends from random fluctuations of the true history. (For purposes 
of the reported investigation, data segmentation into frequency ranges and test 
for time invariance of the sample mean square values in each frequency internal 
was not considered necessary.) 

Amplitude probability analysis may now conveniently be used to evaluate the 
second and third basic characteristics of random data: the presence or absence 

of periodic components, and the normality of the random data. (Power spectral 
density analysis and autocorrelation may additionally be employed for confirma- 
tion of conclusions pertaining to the periodicity reached on the basis of am- 
plitude distribution, or may be used (perhaps less effectively) alone. For the 
BRU simulator, displacement power spectral density analyses were made for all 
random signals at two vibration input levels, as discussed in the next section. 

The results of the probability analysis conducted on the BRU simulator are pre- 
sented in Figures 14 through 23 for the case where external random vibrations 
were applied in the rotor axial direction, and in Figures 24 through 34 for the 
case where external random vibrations were applied in a direction perpendicular 
to the rotor axis . The probability analyses as shown in the above listed figures 
was obtained from the recorded data for which the rms amplitude values had already 
been presented in the series of Figures 77 through 89 and 130 through 142 in 
Reference 1. 

Each of the probability analysis figures is composed of four graphs. Two were 
obtained at the lowest external random vibration input test level of 0.54g (rms) 
and the other two at the highest available test level, which was 1.52g (rms) for 
most plots, but only 1.35g (rms) for a few under transverse external vibration 
input. For each test level, the upper plot depicts the probability density func- 
tion and the lower graph the probability distribution function. The probability 



density curve may be interpreted in the following way: A selected amplitude 

range equal to one two-hundredth (0.02) of half of the total amplitude distribu- 
tion range (corresponding to one window width or bin of the probability analyzer) 
had a probability of occurrence as indicated by the ordinate. The probability 
analysis does not produce a smooth curve because the amount of random data 
analyzed is not infinite. For the case of the BRU data presented, either 328 or 
164 measurements were collected in each of the 400 bins of the analyzer. When 
displayed on an oscilloscope screen (see Figure 35), the 400 probability values 
associated with each of the 400 bins can be individually distinguished. However, 
plotting by an x-y recorder generally produces 'noisy 1 looking probability density 
curves, from which only estimated probabilities can be determined. This is not 
as much of a deficiency as it may at first appear, because the main value of the 
probability density curve may be seen in its visual presentation of symmetry 
about the mean amplitude value. Symmetry about the mean is a necessary but not 
sufficient condition for a 'normal' or Gaussian distribution. The nonsymmetrica 1- 
ly shaped curve is not qualified for simpler statistical treatment, but it still 
may provide some information which immediately relates to the behavior of the 
physical part for which it was recorded. For the data presented in this report, 
moving toward the right-hand side of the probability density curve along the 
abscissa means an increasing signal (voltage) magnitude, which in turn is indica- 
tive of an increasing gap at the location of the capacitance probe in which the 
signal was generated. 

Keeping this physical relationship in mind, examination of the probability density 
curves generated for the rotor-bearing system parts of the BRU simulator indicates 
the following: 

1. Response amplitude distributions are less distorted (more nearly 
normal) at lower external vibration input levels. 

2. Response amplitudes are higher when external vibrations are applied 
- in-line with the direction of motion of the particular part (free 

axis) than perpendicular to it. For example: (a) at the low vibra- 

tion input level, the thrust bearing film thickness variations 
(Figures 22 and 32) are higher when the direction of vibration 
coincides with the rotor axis, than for vibrations applied 



transversely to it. At the high vibration input level, the 
difference in thrust bearing film thickness variations has been 
masked by the fact that in either case the thrust bearing runner 
had already started to move back and forth across the full thrust 
bearing clearance before the highest vibration input level had 
been reached. The gimbal motions, however, are higher for both 
low and high vertical vibration input levels applied in the 
rotor axial direction; (b) when external vibrations were applied 
in the transverse direction, the journal bearing pad to shaft 
film thickness variations as well as the rotor motion relative to 
the casing were larger than for the test case of axially-applied 
external vibrations. Since the maximum amplitude spread of both 
pad flexures (flexure probe 22 for compressor-end pad A and 
flexure probe 21 for turbine-end pad D, Figures 28 and 29 re- 
spectively) increased by more than the increases for their re- 
spective pads plus that of the opposing solidly-mounted pads, it 
must be concluded that all pivots experienced momentary unloading 
and possible pivot-seat separation. (This apparently was not the 
case at any test input level when external vibrations were applied 
axially . ) 

The distortion (skew) of the amplitude response distribution can 
generally be related to the applicable gas film nonlinearity and 
the physical boundaries that inhibit the free motion of the 
vibrating pars. For example: (a) the thrust runner (Figure 22) 

vibrated at low external vibration input levels initially about 
its steady-state operating position, which is approximately 
0.0005 in..away from the compressor-end thrust plate (see Figure 10 
for a schematic of the steady-state thrust bearing operating 
condition). The nonsymmetrical thrust bearing gas film stiffness 
about the operating point permits the rotor to make some very 
large excursions in the direction of the turbine end thrust plate. 
The number of such excursions is small, however, and the result 
is a highly unsymmetrical distribution in the occurrence of am- 
plitudes. At higher vibration input levels, the thrust bearing 
film thickness distribution curve loses all resemblence to the 



Gaussian distribution, and bears close resemblence instead to 
the probability distribution for a sine function. As already 
recognized from the time-amplitude records, this is indicative 
of the thrust runner bouncing back and forth between the two 
stator plates. Figure 22 shows, however, that the probability 
of the rotor being nearer the compressor-end thrust plate is 
higher than being nearer the turbine-end thrust plate; (b) 
another example of highly skewed amplitude distribution curve 
at ^higher vibration levels is given by the film thickness varia- 
tions between the solid-mounted turbine-end journal bearing pad 
and the rotor. Figures 16 and 27 clearly show that the pad is 
vibrating, about a mean position which is approximately 0.G0035 
inch away from the point beyond which no amplitudes occur (the 
rotor surface). The solid-mounted compressor journal bearing 
pad (Figure 26), which was operated with hydrostatic pressuriza- 
tion during the vibration tests, and thus had a higher film 
thickness, does not exhibit similar behavior. For the compressor 
journal pad, the amplitude distribution curve is much more sym- 
metrical, apparently because the pad has more room to vibrate 
(larger steady-state film thickness). Interestingly, its range 
of vibration amplitude is smaller than that of the solid-mounted 
turbine journal pad. However, for the flexibly mounted compressor 
journal pad (which was operated hydrostatically) , the range of 
vibration was larger than that of the hydrodynamically operating 
JElexibly-mounted turbine-end journal bearing pad. _ Since both 
flexures vibrated with nearly equal amplitudes, the 'harder' 
hydrodynamic gas film stiffness of the solid-mounted turbine 
journal pad seems to have provided more 'bounce' to the shaft, 
throwing it back further towards the flex-mounted pad; (c) the 
skew of the amplitude distributions for the two flexures (Figures 
2& and 29) with high vibration input levels in the rotor trans- 
verse direction indicates vibration which may include rotor conical 
motions. This is because at the compressor end, the flexure ' 
vibrates predominately away from the probe while at the turbine 
end most of the vibration occurs closer to the probe. This 
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observation seems to tie in with the previously noted out-of-phase 
vibration of the two flexibly-mounted pads under hydrostatic 
operating conditions, as shown in Figure 7, top photograph. 

The probability distribution function, sometimes also called the cummulative 
probability distribution, is the integral of the density function. The proba- 
bility density analyzer which was used to compute the density function employed 
a bin-by-bin summation process to arrive at the probability distribution. Con- 
sequently, the summation is exact and a smooth curve was produced. Since the 
ordinate of the distribution function is always normalized for 100 percent prob- 
ability, it is very convenient to compare the experimental data curve with a 
superimposed nominal Gaussian distribution function. The closeness of agreement 
between these two curves should give a good visual indication of the normality 
of the data. For easy reference, the 16, 50 and 84 percent points have been 
shown on the ordinate. These correspond to — a, 0, and + a points on the abscissa 
for the Gaussian distribution. The marks on the abscissa which divide the dis- 
tribution into plus and minus three times the standard deviation (+ 3a) do apply, 
of course, only to the Gaussian distribution. For most of the data, the amplitude 

l 

distribution will not equal + 3 a° This was already anticipated from comparison 
of the rms values with the amplitude envelope values as obtained from oscilloscope 
screen photographs, and tabulated for the BRU response data to axially-applied 
external vibrations in Table V. For a number of density and distribution graphs, 
the discrepancy in the abscissa (deviation from + 3a) has been enlarged through 
the plotting process, particularly for very drawn out tail ends and unsymmetrical 
distributions. In those cases the tail ends were cut off and a somewhat magnified 
plot produced in the interest of maintaining a suitable ordinate resolution. 
Nevertheless, comparison between the amplitude values given as plus-minus figures 
in the probability plots, and those obtained from amplitude envelope photographs, 
generally gives agreement within 20 percent. The only significant exception was 
found in the amplitude distribution for the flex-mounted turbine journal bearing 
pad, for which the rms-to-photo comparison had indicated an amplitude spread of 
+ 4a at low external vibration level input (see Table V). The amplitude dis- 
tribution plot indicates a value much closer to what had been obtained for all 
other probes and it is in itself also fairly close to normal (Gaussian). The 
photograph from which the + 4c value, had been obtained would probably have 



contained a rare, excessively large excursion, which by itself was not enough to 
show up on the probability plots. Alternatively, it may not have been included 
in the particular stretch of data that was analyzed for the probability distribu- 
tions. (This incident, which would have been termed 'out-of-range data’ or 
'stray-point' in the analysis of steady-state periodic data, is well within the 
range of possibilities — even if only a very small probability is assigned to 
it — that must be accepted for random data.) 

The probability plots indicate in a few cases that symmetry in the density dis- 
tribution is no guarantee for complete data normality. (See, for example. 

Figure 25.) These plots, however, generally show only slight deviations from 
the Gaussian. Those for which the density function shows significant skew or 
asymmetry cannot be treated by Gaussian statistics. The amplitude analysis 
test results give good indication that for certain types of moving components 
of a film-type bearing system, Gaussian response amplitude distribution cannot 
be expected for any but very small external vibration inputs. 

The following table classifies the various BRU bearing system part responses as 
nearly Gaussian, approximately Gaussian and non-Gaussian. These classifications 
are arbitrary and are included only to serve the purpose of categorizing the 
data. 

Inspection of Table VI clearly indicates a shift towards non-Gaussian distribu- 
tions with the external vibration inputs applied in the transverse direction. 

Such a conclusion does not come as a surprise, however, because as has already 
been observed, transverse vibration inputs elicit much stronger responses from 
the journal bearing pads than do axial vibration inputs. Since the reduction of 
axial vibrations under transverse inputs is not sufficient to make the film 
thickness variation between thrust runner and thrust plates behave in a Gaussian 
distribution, the addition of the non-Gaussian distribution from the journal pads 
and flexures (vibrating more severely with this direction of input) makes the 
total number of non-Gaussian distributions go up. 



- 29 - 




Frequency Analysis — Displacement Power Spectral Density 


Perhaps the most important single descriptive characteristic of stationary random 
data is the power spectral density function (PSD), which defines the frequency 
composition of the data . Under random vibration conditions, the occurrence of 
any particular frequency at which an object moves at any given point in time 
follows a random variation. In this respect it may be compared to the amplitude 
variation, which also follows a random pattern. The significant difference is, 
however, that the distribution of frequencies cannot be estimated or even guessed 
by visual inspection of time-amplitude data only. For constant parameter linear 
physical systems, the output power spectrum is equal to the input power spectrum 
multiplied by the square of the gain factor of the system. For a unity gain 
factor the expected response should therefore equal the frequency distribution 
of the external vibration input. In the case of bandwidth- limited input at con- 
stant acceleration, the system response should show equal accelerations at all 
input frequencies, in effect producing a rectangular acceleration distribution 
over the input signal bandwidth (frequency range). 

For the BRU, the specified external vibration input acceleration is not constant 
over the specified frequency range from 20 to 2000 Hz, but increases logaritmical 
ly from 20 to 100 Hz and then remains constant for the remainder of the spectrum 
(see Figure I). Figure 36 shows the same specification test power spectral 
density converted from acceleration power spectral density to displacement power 
spectral density, (double amplitude displacement). The displacement power spec- 
tral density response data from all bearing system components in the BRU simulate 
should at least resemble Figure 36 if their responses were linear, which essen- 
tially means being resonance- free . 

Information relative to the dynamic characteristics of the system is the primary 
result of the power spectral density analysis. A secondary application includes 
its use for the detection of periodicities . A highly resolved spectrum will 
reveal periodic components as sharp peaks, even when the periodicities are of 
relatively small intensity. Since a sharp peak in the power spectrum of the 
sample data may also represent narrow-band random data (resonances), attempts 
must be made to distinguish these two cases from one another through repetition 
of the power spectral density measurement with finer resolution (narrower -filter 
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bandwidth) . If the measured spectral peak represents a sine wave, the indicated 
bandwidth of the peak will always be equal to the bandwidth of the analyzer 
filter, no matter how narrow the filter. Furthermore, the indicated spectral 
density, will always increase in direct proportion to the reduction in filter 
bandwidth. (This method of detection will clearly not work unless the resolu- 
tion bandwidth of the analyzer filter is smaller than the bandwidth of possible 
narrow-band random data. It is also necessary that the available sample record 
be sufficiently long to make the product of filter bandwidth and analysis time 
(in seconds) much larger than unity. Otherwise, uncertainties in the spectral 
density measurements of the random portions of the data will produce spurious 
peaks which may obscure the peak caused by the sine wave.) 

Examination of the BRU displacement power spectral density plots does not reveal 
any major periodicities in the intermediate test frequency range where most of 
the system resonances occurred. This was ascertained through multiple data 
analysis runs of typical stretches of random data with reduced filter bandwidth 
settings on the real-time analyzer. For example, in Figure 37(c) the same basic 
random data was used for analysis as for Figure 37(a), but the filter bandwidth 
was reduced from 6 Hz in Figure 37(a) to 1.5 Hz in Figure 37(c). The result is 
a much more finely structured analysis in the frequency domain, which, however, 
is limited to the total range extending from 5 to 500 Hz. (Basic instrument 
limitations cause the upper analysis range to shrink from 2000 Hz at 6 Hz filter 
bandwidth to 500 Hz at 1.5 Hz filter bandwidth and finally to 100 Hz at 0.3 Hz 
filter bandwidth.) Further reduction in the filter bandwidth to 0.3 Hz (Figure 
37(d)) does not reveal any major deviations from the essentially flat response 
as already indicated in Figure 37(a). There is, however, a fairly weak peak 
around 60 Hz, which may indicate a sinusoidal component in the BRU system re- 
sponse (possibly caused by line frequency shaker stray input). Alternatively, 
it may only signify a small line frequency voltage superposition in one of the 
measurement'; circuits;.- . The; results : shbwn;'inr.Figuree37'/^etev.obtaihed.; for one ■ 
of the solid-mounted pads in the compressor journal bearing. They seem to be 
typical for system components in BRU simulator, with the possible exception of 
the thrust bearing gimbal and the flexibly-mounted journal bearing pad and 
flexure responses. Narrower filter bandwidth resolution analysis for these 
parts indicates the presence of previously hidden peaks and around 60 and 180 Hz 
(Figures 59, 61 and 75). 



There is one secondary application of the frequency analyses that arises from 
the inherent capability of the real-time analysis instrument: the estimate of 

data stationarity . Data stationarity has already been examined as a prerequisite 
step to amplitude analysis. The stationarity was deduced from the available rms 
plots (Figures 12 and 13, and also from Reference 1). During frequency analysis, 
data stationarity may also be estimated by the following process: With the 
analyzer set for a maximum frequency of 2000 Hz and the sampling of one ensemble, 
the instrument will collect one data point from the random input that falls into 
each of the 500 'bins.' (For the 2000 Hz range setting, each bin is normally 
4 Hz wide.) After having accepted just one point into each of the 4 Hz wide 
bins, the instrument memory is full and hence rejects all further inputs. The 
analyzer is now ready to plot the power spectral density versus frequency curve 
from its memory. The above process may then be repeated from the same random 
data but with increasingly higher numbers of data points collected into each 
analyzer bin. The analyzer is capable of accepting a maximum of 4096 data points 
per bin (4096 ensembles) for which a stretch of random input data 17 minutes long 
would be required. For the BRU simulator, enough test data was taken at con- 
stant vibration input levels for a maximum of 256 ensembles. Power spectral 

| 

density plots for various parts of the BRU simulator were therefore obtained 
for 1, 8, 64, 128 and (if possible) 256 ensembles. These plots show clearly 
.that the composition of the recorded BRU random response data does not notice- 
ably differ for very short stretches of data from those obtained for much longer 
segments. A typical series of power spectral density plots obtained from the 
same data record, but with increasingly higher ensemble rates is presented in 
Figure 37(b) (8 ensembles), Figure 37(a) (64 ensembles), and Figure 41 (256 
ensembles). 

For the BRU simulator, the power spectral density analysis (Figures 38 through 
67) give the following main results: 

- 1. For external vibrations applied in the rotor axial direction, the 
thrust bearing film thickness variation indicates a nearly unity 
gain over the frequency range of approximately -100 to 1000 Hz. 

(See Figures 54 and 55.) This represents the highest 'base line' 
response measured (excluding the resonance peaks) in either 
direction of applied vibrations and for any recorded response. 
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2. Other bearing-system components with high response-to- input 
ratios (which actually amounts to low attenuation, because the 
gain factor is less than one) for the frequency range 100 to 1000 
Hz are the thrust bearing gimbal motions (vibration applied in 
the rotor axial direction) and the film thickness variation be- 
tween the solid-mounted pads and the rotor (with vibrations 
applied in the rotor transverse direction). 

3. In the frequency range between 20 and approximately 60 Hz, the 
response of all bearing system components is either at almost 
constant displacement power spectral density level or exhibits 
slightly increasing levels with frequency. Considerable attenua- 
tion of the response vibrations relative to the shake table input 
vibration in this frequency region must be noted. However, the 
vibration response level extends into the low frequency region 
(5-20 Hz) for which no nominal excitation was present. The 
level of vibration response in the region up to 60 Hz shows a 
definite, positive, statistical correlation to the input level. 

4. Significant power spectral density response peaks (which exceed 
the input power spectral density for any of the parts investi- 
gated) occurred between 100 and 800 Hz. Except for the response 
peaks at 600 Hz, the general distribution of resonances matches 
previously obtained results for sinusoidal vibration input 
scanning. The very narrow resonance peak at 600 Hz obviously 
corresponds to the running frequency of the BRU simulator and 
most probably represents the sinusoidal vibration due to resid- 
ual rotor unbalance. (Sinusoidal vibration testing did not reveal 
any resonance at 600 Hz beyond the once-per-revolution rotor un- 
balance response.) From evaluations of the sinuosidal vibration 
response results obtained from the BRU simulator, it may also be 
concluded that the response peak amplitudes at 800 Hz are too 
small in absolute values to be of concern. 

The available experimental data does not yield positive identifi- 
cation of the gas-film resonance frequencies associated with the 


5 . 


bearing ' system . There are at least five such frequencies, some 
of which are located very close together. The broad based peaks 
observed in the power spectral density analyses are probably due 
to these resonances. (The -five resonance frequencies originate 
from the thrust bearing gas film and from the journal bearing 
gas film in the in-line direction with the flexibly-mounted pad 
(direction of minimum gas film stiffness) and at 90 degrees to 
it (direction of maximum journal bearing stiffness). Since the 
compressor journal was operated hydrostatically and the turbine 
journal was operated hydrodynamically during the random vibration 
tests, their stiffnesses are different. Consequently, each condi- 
tion will prpduce : a ; set. of ; resonances :at slightly different fre- 
quencies.) 
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ANALYSIS OF VIBRATION RESPONSE WITH SINUSOIDAL INPUTS 


A considerable amount of sinusoidal vibration response data has been recorded 
for the BRU simulator (Ref. 1). Since complete sets of data were acquired 
under various BRU operating conditions, comparison between sets may be used 
to . draw important conclusions as to the behavior of the BRU 1 s internal compohe.nts. 
This section describes the expected effects of actual test operating condi- 
tions upon this behavior and formulates preliminary explanations of the ob- 
served results. 


Vibration data -is available for the following BRU operating conditions: 
Direction of 


Externally 
Applied 
Vibration 
Level (g) 

Applied Vibration 

Operating Mode of BRU 

Bearings 


Along 

Rotor 

Axis 

Across 

Rotor 

Axis 

Compressor 

Journal 

Bearing 

Turbine 

Journal 

Bearing 

Thrust 

Bearing 

Figure Numbers 
in Reference 1 

.25 

X 


Hydro- 

static 

Hydro- 

static 

Hydro-.. 

static 

, 42-52 

.12 


X 

Hydro- 

static 

Hydro- 

static 

Hydro- 

static 

105-116 

.12. 

X 


Hydro- 

dynamic 

Hydro- 

dynamic 

Hydro- 

dynamic 

53-71 

.12 


X 

Hydro- 

static 

Hydro- 

dynamic 

Hydro- 

dynamic 

117-129 


The two different bearing operating modes are important to the interpretation 
of the vibration response data in that the bearing gas; film spring stiffnesses 
are different for each mode. For hydrostatic bearing operation, the journal- 
bearing gas film thickness is usually larger than that for hydrodynamic op- 
eration. Therefore, the bearing usually has lower stiffness in the hydro- 
static operating condition. This may be expected to produce rotor bearing 
system critical speeds at lower frequencies for hydrostatic operation than 
for hydrodynamic operation. .. For gas-lubricated thrust bearings, the opposite 
may be true. Since total thrust bearing clearance is fixed (in contrast to 



the spring-loaded journal bearings), the relatively high gas supply pressure 
for hydrostatic operation tends to produce a fairly high thrust bearing gas 
film stiffness. The manner in which the thrust bearing was operated hydro- 
dynamically during the vibration tests makes its gas film stiffness somewhat 
different than that expected during normal operation: the rotor was supported 
against gravity by the differential upward pressure in the simulator housing, 
to the extent that the net runner loading was against the compressor (upper) 
thrust plate (see Reference 1 , page 24-25). The gas film thickness between 
the compressor-side thrust plate and the runner was approximately 0.0004 to 
0.0005 inch, and thus* this side contributed the majority of the axial stiff- 
ness. The stiffness of the turbine side (lower) thrust bearing gas film 
was probably considerably lower than that of the upper film, because of the 
much higher film on the lower side. Moreover, the end of the rotor was ex- 
posed to the large-volume, low-pressure gas cushion in the circulator cavity, 
an element which contributes very little stiffness. 

A partial structural failure of the compressor-end journal bearing early in 
the test series made it necessary to operate that bearing hydrostatically 
for the majority of the tests. For the case of laterally-applied vibration 
(Figures 117-129 in Reference 1), the combined operation of the hydrostatically- 
pressurized compressor-end journal bearing with £he...hydr!bdynamically-operating 
tufbine^end ' jodtnal bearirig^and thruSt bearing makes: it rather- likely, that ac- 
tual resonance frequencies of the- BRU . simulator- xompbnentsrwil lube different 
from- those which would be expected during hydrodynamic operation of all bearings. 

While the actual resonance frequencies (and mode shapes) are affected by the 
necessity of operating one journal bearing hydrostatically, the rotor may 
still be expected to exhibit some evidence of resonance (probably of low amp- 
litude since the rotor is well balanced) at two distinct speeds. 

The mode shapes at these speeds will be essentially either cylindrical or 
conical in nature. The order of occurrence of these modes will, of course, 
depend upon relationships between rotor mass, bearing span, polar amount of 
inertia, and moment of inertia about the rotor center of gravity. For the 



BRU rotor with fully hydrodynamic bearings, the first mode has been calculated 
as conical at 5,270 rpm (88 Hz) and the second mode as essentially translatory 
at 5,680 rpm (94 Hz). The frequency ratio between the conical and the trans- 
latory modes is thus 1.07 for the fully-hydrodynamic case (see Reference 2, 
page 97). 

The tilting-pad gas-lubricated journal bearing, in the three-pad configuration 
found in the BRU, has distinctly different radial spring stiffnesses in the 
direction of the pivot axis of the flexure-supported pad (referred to as the 
y-direction) compared to that in the direction perpendicular to it (referred 
to as the x-direction) . Ratios between 2.5 and 6 for the stiffness in the 
x-direction to that of the y-direction seem to be typical for this type of 
bearing at commonly-operated gas film thicknesses. (Reference 3, page 4-14, 
lists calculated spring stiffness ratios of 2.8 for a number of operating con- 
ditions. Calculations performed at MTI for a journal bearing with dimensions 
and test operating conditions identical to those in the BRU show ratios from 
3.4 at 0.00025 inch minimum film thickness to 6 at 0.0005 inch film thickness. 
Other, unpublished test results obtained by the senior author for a very 
similar gas journal bearing indicated a similar bearing stiffness component 
ratio of 2.5) Since the frequencies of both modes of rigid body vibration 
are proportional to the square root of the bearing stiffness, the frequency 
ratio between vibration in the x- and y-directions reduces to values between 
1.58 and 2.45 for either mode. A typical mean value might be 2.2. 

The foregoing discussion has identified the existence of four lateral vibra- 
tion natural frequencies for the nominal BRU rotor bearing system: conical 

mode y-direction, translatory mode y-direction, conical mode x-direction, and 
translatory mode x-direction. The calculated frequency ratios are for fully- 
hydrodynamic BRU operation: conical y to translatory y, 1:1.07 (88:94 Hz); 

conical y to conical x, 1:1.58 to 1:2.45 (88:139 to 88:216 Hz); and conical 
x to translatory x, 1:1.07 (139:145 to 216:232 Hz). 

t 

Examination of the recorded lateral rotor motion vibration response, with 
external vibration applied in the rotor axial direction and full hydrodynamic 


bearing operation, indicates a resonance pattern resembling the one described 
above. (See Figures 55 through 61 in Reference 1 for lateral rotor motions 
and associated pad-to-shaft film thickness variations.) The lowest of the 
four resonance frequencies occurs at approximately 100 Hz, with the three 
next lowest occurring at approximately 130, 150 and 185 Hz. 

The fact that the lowest resonance frequency is observed at 100 Hz, instead 
of at 88 Hz as predicted, may well be due to the lower journal bearing oper- 
ating film thickness (and thus higher stiffness) in the N 2 atmosphere during 
the vibration tests as compared to the helium-xenon atmosphere for which the • 
resonance frequencies were predicted in Reference 2. The observed ratio of 
approximately 1.3 for conical to translatory resonances indicates a journal 
bearing stiffness ratio of approximately 1.7 for the x and y directions. 

This value is well within the values of 1.58 to 2.45 observed for similar 
bearing-rotor systems as described above. 

The weak resonance peak at approximately 70 Hz in the cited reference figures 
is apparently caused by the thrust bearing resonance. This resonance, as 
exhibited in Figure 62 of Reference 1, is of a somewhat unusual type: a 

fairly gradual buildup of rotor amplitudes (in the vertical direction) fol- 
lowed by a very steep and sudden decline; ■ Ihe same behavior is observed for 
the thrust bearing gimbal (Figure' 63, Reference 1). Since a distinct phase 
shift of nearly 180° has been found to accompany the drop in thrust runner- 
to-stator gas film variations, it is concluded that the thrust bearing re- 
sonance does occur at approximately 70 Hz.- (Figures .62 and 63 in Reference 1 
show this resonance at 60 Hz. Careful reexamination of the tape-recorded data 
indicates, however, that the resonance actually occurred between 69 and 70 Hz. 
The deviation in these two figures in Reference 1, appears to be the result 
of nonlinearities in the x-y plotter, which was calibrated at the end-points 
of the horizontal scale for these particular plots.) 

Since appreciable gimbal motion was found to be associated with the thrust 
runner- to-stator resonance, is apparent that rotor wobble was present. This 
offers at least a partial small peak just above 70 Hz in the plots of lateral 
rotor response. , 



The interpretation of the resonance pattern of rotor motion due to lateral 
external vibration inputs is considerably more difficult. The thrust bearing 
resonance was not expected to have shifted appreciably in frequency. This 
seems to be borne out by Figure 128 in Reference 1, which indicates the major 
peak to lie between 70 and 80 Hz. However, the bearing response now also 
indicates several smaller peaks, which may be the result of lateral rotor 
motions about the rotor center of gravity, induced by the nonsymmetrical : 
mounting of the BRU simulator the the shaker. The lateral rotor vibration 
response pattern is less structured than that observed for the case of axially 
applied external vibrations. Strong resonances now occur just below and just 
above 100 Hz, and a definite downward shift; in frequency is noted. This could 
well be the result of reduced stiffness in the compressor journal bearing, 
which was operated hydrostatically during these tests. 

A distinct resonance peak which occurs quite persistently between 300 and 400 
Hz in most of the response plots obtained under lateral external vibration 
inputs, is believed to be caused by a structural resonance which either ori- 
ginates in the support fixture or in the interaction between the fixture and 
some structural component of the BRU simulator. 


CONCLUSIONS 


Analysis of Response of Linear System to Externally-Imposed Random Vibration 

The extended analysis of the linear three-mass axial response model permits the 
calculation of thrust bearing surface contact probabilities for given (or 
assumed) mean thrust bearing gas-film thicknesses. These calculations are based, 
however, upon the assumption of a normal (or Gaussian) distribution of thrust 
runner- to- thrust plate gas-film thickness variations. 

Experimental Test Evaluation* 

Experimental analysis of the BRU simulator response data in the frequency domain 
(displacement power spectral density) has indicated a considerable degree of 
similarity between the rotor-bearing system responses to externally-applied 
sinusoidal and random vibration inputs. This seems to indicate the absence of 
vibrating subsystems having different resonance frequencies which, if excited 
simultaneously, might tend to reinforce each other to produce severe vibrations 
in the random vibration environment . 

Comparison of the vibration response and input levels (gain factor) indicates 
considerable attenuation for most rotor-bearing system components. Least atten- 
uation was noted for the rotor-thrust bearing vibrations with external excita- 
tions in the rotor-axial direction. 

Amplitude analysis (probability distribution) of the BRU simulator rotor-bearing 
system responses has yielded a very clear picture of the physical behavior of the 
BRU rotor-bearing system components at different levels and directions of 
externally-applied vibration inputs . The experimental probability analysis has 
further shown that the response of most of the gas film supported rotor system 


*The conclusions presented in this report are based upon an analysis and experi- 
mental evaluation of vibration data presented and tentatively discussed in 
Reference 1. They should be considered in conjunction with the conclusions 
presented in Reference 1. 
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components tends to shift toward non-Gaussian distributions as the external 
vibration input level is increased. This appears to be caused by the increasing- 
ly nonlinear spring stiffness characteristics of the gas films at decreasing 
thicknesses. The thrust bearing film thickness variations under rotor axial 
vibration inputs (Figure 22) is a particularly vivid example of this behavior, 
and underscores the need for a practical design tool to predict this phenomenon. 
Further development of the theoretical probability analysis, to include the non- 
linear nature of the gas films, should significantly improve the prediction of 
the liklihood of bearing surface contacts under high levels of externally-imposed 

random vibration. Such development should be included in future investigations 

V 

in this area . 

A detailed examination of the BRU simulator response data obtained under sinusoi- 
dal vibration inputs indicates that the relatively large number of vibrating sub- 
systems (rotor mass - gas bearing film combinations) all have the potential to 
produce severe resonances. For this reason, efforts to suppress the general 
vibration response of the BRU bearing system will have to be broadly based to 
treat the worst condition. 

Finally* it should be noted that the thrust bearing has not had a complete 
evaluation under random vibration due to the nonsimilarity of dynamic gas forces 
on the impellers and cavities in the simulator, compared to the complete BRU* 

In the simulator the rotor was partially supported by gas forces in the cavity, 
which led to a rather low and unrepresentative stiffness of the thrust bearing 
gas film. The results obtained to date for the thrust bearing have, however, 
yielded much of the information needed to assess its dynamic behavior under 
random vibration, including the conclusion that some method of axial vibration 
attenuation must be found. Experimental evaluation of such methods in future 
programs should include close simulation of the actual thrust bearing properties. 
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MATHEMATICAL SYMBOLS 

damping coefficient, lb-sec/in. 
mean square value of ( ) 
frequency, Hz 

mean gas-film thickness, in. 
transfer function 
spring stiffness, lb/in. 
mass, lb-sec /in. 
probability density function 

probability of thrust bearing surface contact 

power spectral density 

input power spectral density 

coordinates defined in Figure 77 

relative coordinates defined in Equation (A-4) 

defined in Equation (A-20) 

mean value of y^ 

standard deviation of y^ 

root-mean-square value of y^ 

frequency, rad/sec 
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APPENDIX A 

PROBABILITY ANALYSIS OF DYNAMIC RESPONSE OF 
BRAYTON ROTATING UNIT TO AXIAL RANDOM EXCITATIONS 

The dynamic response of the BRU thrust bearing due to axial random excitation 
was analyzed in Appendix C of Reference 1. For easy reference,, the analysis is 
recapitualted herein,, together with the thrust bearing gas film thickness varia- 
tion. probability, analysis . 

The Brayton rotating unit was modeled as a three-mass system, as shown in Figure 
77. A sketch of the BRU is shown in Figure 78, in which the various masses and 
springs are identified with the model shown in Figure 77. The base of the three- 
mass system is the shake table, by which random excitations are transmitted to 
the BRU. The quantities k^ and c-^ are, respectively, the stiffness and damping 
of the connectors which provide the joint between the BRU and the shake table; 
m^ is the mass of the BRU housing structure; k^ and c^ are the thrust bearing gas- 
film. stiffness: and damping, respectively; and finally, m^ is the rotor mass. The 
equations; of motion of the three-mass system shown in Figure 77 are: 

m- Xj_ + c ] ,(x 1 -x b ) + c 2 (x 1 -x 2 ) + k 1 (x 1 -x b ) + k^x^x^ = 0 (A-l) 

m 2 X 2 + c 2^2" X 1^ + c 3^2 _X 3^ + k 2^ x 2 _x l^ + k 3 ^ X 2 _X 3^ = ° ( A_2 ) 

m x’ 3 ..+ c 3 (x 3 "X 2 ) + k 3 (x 3 ~x 2 ) = G (A-3) 

Introducing, relative coordinates between mass points, 

Y i. = X 1 ' X b 

y 2 = X 2' " X 1 (A-^) 

y 3 - x 3 - x 2 



Equations (A-l) and (A-3) become 


m l y l + C l y l ' C 2 y 2 + k l y l " k 2 y 2 = “Vb 


(A-5) 


m 2 (y i + y 2 } + C 2 y 2 ' c 3 y 3 + k 2 y 2 ' k 3 y 3 = ' m 2 X b 


m 3 (y l + y 2 + y 3 ) + C 3 y 3 + k 3 y 3 = ~ m 3*b 


(A-6) 
(A- 7) 


Dividing Equations (A-5), A-6) and (A-7) by m^ , and respectively, 


C i C k , k o 

, 1 • 2 • , 1 2 

y, + r - y.i ■ “ y? + ~ yi ■ ~ y 9 = . 

i m^ 1 2 m^ L m^ z b 


(A-8) 


•• , . C 2 . C 3 • , k 2 k 3 

y l + y 2 + ^ y 2 " ^ y 3 + y 2 ‘ SJ y 3 = " X b 


(A-9) 


c o k * 

y, + y 9 + + — y. + — = -x 

l 2 J 3 m 3 3 m 3 3 ° 


(A- 10) 


Subtracting Equation (A-8) from Equation (A-9), 


r c_ c_“i 

c„ k. 

r k o 


2,2 

• 3 * 1 

2 . 

2 

h 

y 0 - — y 0 - — y. + 

— + — 

L m l m 2_ 

J 2 m 2 ^3 m^ J 1 

l m l 

m 2 J 


y 3 = 0 (A- 11) 


Subtracting Equation (A-9) from Equation (A-10), 

? + r^ + ^i 

y 3 m 2 y 2 [_m 2 m 3 J y 3 m 2 y 2 [_m 2 mj 


y 3 = ° 


(A- 12) 


Equations (A-8), (A- 11) and (A-12) can be written in the matrix form as follows: 





(A-14) 


To obtain the complex frequency response, assume 


x b = e 


H. (to) 


y 2 = H v («) ‘ e 

y i 


(A-15) 


H (to) 


and substitute into Equation (A-13) : 


|H (to) 


fH (to) 


H (to) 
y l 


- u [I] H (to) + ito[C] H (to) + [KJ H (to) 

y 2 y 2 y 2 


= 0) 0 (A-16) 


H (to) 
y 3 


H (w) 


H (to) 
y 3 


1 0 0 


0 1 0 


( A-17) 


0 0 1 


^2 + ^2 
m l m 2 


(A-18) 


f3 + ^3 
m 2 m 3 
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Define a matrix Z: 

[Z] = - [I] + ii [C] + ~ [K] 

01 l 

0) 

Equation (A-16) can be written as 


H 

(«)“ 


Y 

y l 



H 

(w) 

2 

= CO 

0 

y 2 



H 

(ui) 


0 

l y 3 





Thus , 


H 

( to) 


1 ' 

y l 



H 

(«) 

= CO 2 [Z]" 1 

0 

CM 

>> 



H 

(10) 


0 

. y 3 





(A- 19) 


(A-20) 


(A-21) 
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The complex frequency responses, H (u>) , H (u>) and H (w) can be readily 


calculated from Equation (A-21). The complex frequency response functions are 
sometimes called transfer functions. They give the response at one point due to 
a unit input excitation at some other point. 


It was shown, for example, in References 4 and 5 that the input and output power 
spectral densities of a random process are also related by the transfer functions 
in the following manner. 


S («o) 
y i 


I h M 2 

y i 


S Q («) 


(i=l, 2, 3) 


(A-22) 


where 


S («) 
y i 


s o (to) 

H (w) 
y i 


PSD of the system output y^ (i=l, 2, 3). 

PSD of the system input. 

magnitude of transfer function as calculated from (A-21) . 


Of particular interest is S (cn) which is the power spectral density of the thrust 
bearing gas-film thickness variation. The mean square value can be obtained by 
integrating the power spectral density throughout the whole spectrum. 


cn„ 


y 3 


- s 


S-. (co) du> 


(A-23) 


where \jf^ = mean square value of y^ 


co. 


= lower and upper limits of the spectrum 


The mean square value is related to the variance and the mean value p by the 
following: 



(A-24) 



Since- y 2 is the thrust hearing gas-film thickness variation, its mean value is 
zero- for steady-state operation. 

(i =- 0 (A-25) 

y 3 


Therefore?, knowing the power, spectral density-of y^ we calculate from Equation 

O 

(A-23) to (A-25) , the: standard deviation or. the variance Cy^ . 


2 

CTy- 

y 3 


= t 


yj 




CO, 


Sy 3 (co) dco 


(A-26) 


The random excitation output of the shake table has a normal or Gaussian dis- 
tribution. The. shake table is rigidly connected to the BRU base which is shown 
in. Figure 77 . Thus:,, the: input axia l random exc itation of the BRU is the same as 
the- shake: table, output and has a Gaussian distribution. 

Now;,, the model, of the BRU. rotor-bearing system, as shown in Figure 77,. is a 
linear: constant-parameter system:. For such a system, a random input with 
Gaussian distribution would yield a random response with the same distribution 
(References: 6 and 7) . 

Thus ,, the. Gaussian probability density function, of y^ can be expressed as 


Ply 


<?y 3 V^fr 


* exp 


(A-27) 


2 a 


*3 


Let: h be the. mean, gas:- film thickness . Then, whenever y^ is greater than or 
equal to h, there is: solid contact or rubbing between bearing surfaces. The 
probability that y^ > h or the probability of bearing surface contact, is given 
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K 


where y 3 = y 3 /a y ^ 

TT = h/ff 

y 3 

The values of integral in (A-28) have been tabulated versus h, for example on 
page 387 of Reference 7. 



B-l 


APPENDIX B 

COMPUTER PROGRAM FOR PREDICTING THE DYNAMIC RESPONSE 
OF THE BRAYTON ROTATING UNIT TO AXIAL RANDOM EXCITATIONS 
AND PROBABILITY OF THRUST BEARING SURFACE CONTACT 


This Appendix presents a complete listing of the computer program MTI-492, 
which calculates the axial response of the Brayton Rotating Unit under random 
excitation and the probability of thrust bearing surface contact. Also in- 
cluded in this Appendix are instructions for preparing input data cards, sample 
input and sample computer output. 
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PROGRAM INPUT CARDS 
Card 1 FORMAT (15) 

This card reads in NCASE which is the number of cases. 

Card 2 FORMAT (80H) 

This card contains the descriptive identification for the case being computed. 

Card 3 FORMAT (8E10.3) 

This card reads in FAC. 

FAC = numerical factor to be multiplied with the power spectral density 
specifications to permit calculations at various levels without 
requiring the entire specification to be modified and resubmitted as 
input. 

Card 4 FORMAT (8E10.3) 

This card contains ENGM(l), ENGC(l) and ENGK(l). 

ENGM(l) = m^ in pounds 

ENGC(l) = c^ in lb sec/in. 

ENGK(l) = ^ in lb/in. 

The quantities m^, c^ and are shown in Figure 77. 

Card 5 FORMAT (8E10.3) 

This card contains ENGM(2), ENGC (2) and ENGK(2). 

ENGM(2) = m2 in pounds 

ENGC (2) - c 2 in lb sec/in. 

ENGK (2 ) = k 2 in lb/in. 

The quantities m 2 , c 2 and k 2 are shown in Figure 77. 

Card 6 FORMAT (8E10.3) 

This card contains ENGM(3) , ENGC(3) and ENGK(3) . 

ENGM(3) = m^ in pounds 

ENGC (3) - c 3 in lb sec /in. " " 
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ENGK(3) = k 3 in lb/ in. 

The quantities m^, and are shown in Figure 77. 

Card 7 FORMAT (515) 

This card contains NRANG, SPEC, NFILMT. 

NRANG = number of frequency ranges. Each frequency range will be defined in 

Card 8 by an initial frequency, a final frequency and an incremental 

frequency. 

SPEC = 1: PSD of Figure 1 will be read in. 

2 

2: a white noise displacement PSD of 1 in /Hz will be read in. 

2 

3: a white noise acceleration PSD of 0.01 g /Hz will be read in. 

NOTE: The PSD values read in will be multiplied by FAC specified 

in Card 3. 

NFILMT = Number of film thicknesses (up to a maximum of 10). 

Card 8 FORMAT (8E10.3) 

This set of cards specifies the frequencies at which the response PSD will be 
calculated. The number of cards in this set is equal to NRANG of Card 7. Each 
card containing, WINIT, WFIN and WDELT represents a frequency range. 

WINIT = initial frequency in Hz. 

WFIN = final frequency in Hz. 

WDELT = incremental frequency in Hz. 

Card 9 FORMAT (8E10.3) 

This card contains FLMTH(I), which is the film thicknesses in inches. A total 
of I = NFILMT values is to be read in. 

NOTE: Card 2 through Card 8 represent one case. The whole sequence can be 

repeated as many times as the number of cases specified in Card 1. 
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MECHANICAL TECHNOLOGY INCORPORATED 

FEB173 

1 

968 ALBANY-SHAKER ROAD 

FEB173 

2 

LATHAM, NEW YORK 12110 

FEB173 

3 

518-785-2211 

FEB173 

4 


FEB173 

5 

AXIAL RESPONSE OF TURBOMACHINEPY TO 

FEB173 

6 

STATIONARY RANDOM EXCITATION - OCTOBER* 1972 

FEB173 

7 


FEB173 

8 

DIMENSION TITLE (101 ,M(3,3),C(3,3),K(3,3> 

VIBRAN 

3 

CIMENSION WNA T ( 3 ) *ZETA(3) *ENGM (3) ,ENGC(3) *ENGK(3> 

VIBRAN 

4 

DIMENSION WIN IT (5) .WDELT (5) ,wFIN(5) ,E I N I T ( 5 ) « FOFLT ( 5 ) * FEIN (5) 

VIBRAN 

5 

DIMENSION PsD(3) .PSOY (31 

VIBRAN 

6 

DIMENSION YINj (3) *tOty (3) 

VIBRAN 

1 

DIMENSION RMS ( 5) , Z ALP m A (501 ,PROB(50) 

AU3072 

i 

COMMON/FT/ ELMTh(IO) *nejlmt 

AU 3 O 72 

2 

REAL M.K*ML) 

VIBRAN 

8 

integer spec 

VIBRAN 

9 

C AT A ( ZALPr A ( I ) * 1 = 1 » A 2 ) / , 0 » . 1 » .2* . 3* .A * # 5* .6* .7* .8 * ,9* 1 . * 1 • 1 » 1 #2 * 

AU3072 

3 

1 1.3,1. A, 1.5,1 .6, 1.7,1 .8,1 .9, 2. 0*2.1 ,2. 2*2. 3*2. 4, 2. 5. 2. 6*2. 7*2.8, 

AU3072 

4 

2 2. 9, 3. 0,3.1 *3. 2 ,3. 3, 3. 4, 3. 5, 3. 6, 3. 7, 3. 8, 3. 9, 4. 0,4. 1/ 

AU3072 

5 

DATA (pROB(I) ,1 = 1,42)/ . 5, . 4602, . 4207, . 3821 ,.3446,. 3085,. 2/43, 

AU3072 

6 

1 .2420* *21 19, . 1841, .1587, .1357, .1151, .0968, . 0808 ,. 0668, . 0548, 

AU3072 

7 

2 .0446, .0359, .0287, .0228, .01 79, . 0 1 39 , .0 1 07, . 00820 ,. 00621 ,400466, 

AU3072 

8 

3 .00347, .00256, .00187, .001 30, .00 100,. 00 070*. 00 050*. 00034, 400024, 

AU3072 

9 

4 .00017, .0001 1 , .00006, .00002, .00 , .00/ 

AU3072 

10 

N ALPHA=42 , 

AU3072 

11 

READ ( 5, 1 15 INCASE 

VIBRAN 

10 

KCASE=0 

VIBRAN 

11 

5 CALL REAOIT (TITLE, ENGM,ENGC,ENGK, NRANG, SPEC, W1NIT,WFIN,WDLLT,EAC) 

VIBRAN 

12 

CALL WRIT IT ( TITLE, ENGM*ENGC,ENGK,NRANG, SPEC, WINIT,WF I N.WDtLT »EAC> 

VIBRAN 

13 

CALL MAT RX2 (M,C,K,ENGM,ENGCENGK) 

VIBRAN 

14 

WRITE (6* 195) 

VIBRAN 

15 

DO 2 1=1*3 

VIBRAN 

16 

TOTY < I ) =0. 

VIBRAN 

17 

2 WRITE(6*1S0) (M ( I , J> » J= 1 ,31 , (C ( I , J • , J=l , 3) , (K ( I » J) , J= 1 , 31 

VIBRAN 

18 

TOTS=0. 

VIBRAN 

19 

TOTP=0. 

VIBRAN 

20 

DO 420 1 1 1 = 1 ,NRANG 

VIBRAN 

21 

WRITE <6’90) 

VIBRAN 

22 

WRITE <6,1001 I I I, WI NIT (III) ,WFIN(III» ,WDELT(III) 

VIBRAN 

23 

WRITE (6,110) 

VIBRAN 

24 

W=WINIT < 1 1 1 ) 

VIBRAN 

25 

DEL=WDELT(III) 

VIBRAN 

26 

W=w-DEL 

VIBRAN 

27 

DO 650 3=1,3 

VIBRAN 

28 

650 YINT(J)=0. 

VIBRAN 

29 

PSINT=0 • 

VIBRAN 

30 

SCAL I N=0 . 

VIBRAN 

31 

DO 400 11=1 *1000 

VIBRAN 

32 

W=W+DEL 

VIBRAN 

33 

XW=W 

VIBRAN 

34 

IF (W.GT.WEIN ( I I I ) ) GO TO 4Q0 

VIBRAN 

35 

ALL CALCULATIONS ARE DONE IN RADIANS/SEC 

VIBRAN 

36 

W=W*2« *3. 14159 

VIBRAN 

37 

CALL CALC (M,C,K,W, PSD) 

VIBRAN 

38 

convert all frequencies in to hz before calling for specs 

VIBRAN 

39 

F=W/(2.*3.14159) 

VIBRAN 

40 

Go TO (585,600,610,600,600 >*S PEC 

VIBRAN 

41 

585 CALL SPEC! (F, SCALE) 

VIBRAN 

42 

GO TO 645 

VIBRAN 

43 

600 CALL SPEC2(F, SCALE) 

VIBRAN 

44 

GO TO 645 

VIBRAN 

45 

610 CALL SPEC3 (E, SCALE) 

VIBRAN 

4b 
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&4S 

CONTINUE 


vibran 

47 


scal£=scale*fac 


VIBRAN 

48 


I F ( w » NE • 0 • ) PSINCH=SCAlE*(386.**2/W**4> 


VIBRAN 

49 


IF(W.EQ.O.) PSINCH=SCALE*(386.**2/(1.E-15>**4) 


VIBRAN 

50 


IFdl.EO.l.OR.XW.GE. (WFIN(III)-.Ol) ) GO TO 200 


VIBRAN . 

51 


PSINT=PSINT*PSINCH 


VIBRAN 

52 


scalin=scalin*scale 


VIBRAN 

53 


GO TO 220 


VIBRAN 

54 

200 

PSINT=PSINT*.5*PSINCH 


VIBRAN 

55 


SCALIN=SCALIN+. 5 *SCALE 


VIBRAN 

56 


if(xw.lt. (wfin(iii)-.oi > > go to 220 


VIBRAN 

57 


PSINT = PSIM*W0ELT (III) 


vibran 

58 


SCALlN=SCALIN*WOELT (III) 


vibRan 

59 

220 

CONTINUE 


VIBRAN 

60 


CO 330 J=l»3 


VIBRAN 

61 


PSOY ( J)=SCAlE*PSD( J)*386.**2 


VIBRAN 

62 


IFdl.EO.l.OR.XW.GE. (W^IN(III)-. 01) >GO TO 350 


vibran 

63 


YINt ( J> =YINT <J> ♦PSOY ( J) 


VIBRAN 

64 


C-0 TO 330 


VIBRAN 

65 

350 

YINT ( J) =VINT ( J) ♦. 5 1, PSOY ( J) 


VIBRAN 

66 


IFUW.GE. (WFIN(III)-.Ol ) ) YINT(J)=YINT(J)*WDELT(III) 


VIBRAN 

67 

330 

CONTINUE 


vibran 

68 


CONVERT ENGINEERING FREQUENCIES SACK INTO HZ 


VIBRAN 

69 


W=W/ (2. *3. 14159) 


VIBRAN 

70 


WRITE ( 6 , 120) W, (PSOY(J) ,J=1,3) .PS INCH, SCALE 


VIBRAN 

71 

400 

IF (II. EO. 1000) WRITE ( 6 , 160) < Y I NT < J > , J= 1 * 3) ,PSINT , SCALIN 


VIBRAN 

72 


DO 30 J=1 ,3 


VIBRAN 

73 

30 

TOTY(J)=TCTY(J)+YINT(J) 


vibran 

74 


TOTP=TOTP+PSINT 


vibran 

75 


tots=tots+sca l i n 


vibran 

76 


WRITE (6, 170) (TOTY(J) , J=l,3> ,T0TP,T0TS 


vibran 

77 


CO 40 J=W3 


AU3072 

12 

40 

PMs(J)=sQPT(TOTY(J) ) 


AU3072 

13 


RMS(4)=SQRt (TOTP) 


AU3072 

14 


RMS(5)=SQRT (TOTS) 


AU3072 

15 


WRITE(6,180> (RMS ( J) , J=1 ,5) 


AU3072 

16 


IF(NFILMT.EQ.O) GO TO 420 


AU3072 

17 


IF(III.NE.NRANG) GO TO 420 


0C1872 

i 


WRI TE (6 , 125 ) 


AU3072 

18 


CO 50 J=1 «NF ILMf 


AU3172 

1 


ZAL=FLMth ( j) /RM S ( 3 ) 


AU3172 

2 


CALL tLu<2AL,RPR0R»ZALPHA,PR0B,NALPhA) 


Au3172 

3 

50 

WRItE ( 6,130) FLMjM J> fRPROe 


AU3172 ■ 

4 

420 

continue 


VIBRAN 

78 


KCASE=KCASE >1 


VIBRAN 

79 


IF (KCASE-NCASE)5,425,42S 


VIBRAN 

80 

425 

CALL EXIT 


VIBRAN 

8 l 

90 

FORMAT ( 1 H 1 ,//25X,80 (1H*> //) 


vibran 

82 

100 

FORMAT (25X,22HFREQUENCY RANGE NUM8ER*l5/ 


VIBRAN 

83 

1 

25x,21HCOVERS THE RANGE FROM ,F9 . 2 , 3X , 2HTO , F9 . 2/ 


VIBRAN 

84 

c 

! 25X, 19HIN INCREMENTS OF ,F7 .2/) 


vibran 

85 

110 

FORMAT <25X,80H FREQUENCY PSD OF Y1 PSD OF Y2 

PSD OF Y3 

VIBRAN 

8 b 

1 

psd Input psd input / 


VIBRAN 

87 

2 

: 25X,80H (HZ) (In*In/HZ) (In*In/HZ) 

(lN*lN/HZ) 

VIBRAN 

88 

3 

I ( In* In/HZ ) (G«G/HZ) /) 


VIBRAN 

89 

115 

F ORM AT ( 15 ) 


VIBRAN 

90 

120 

FORMAT (25x,F 12.3,4X,5 (El 0.3, 3 x)» 


VIBRAN 

91 

160 

FORMAT (/23x,16 H I NtEGRAtEO VALuE,2x,5 CE10.3,3X) > 


Au3072 

23 

170 

FORMAT (23X.16HT0TAL INT. V ALUE , 2X , 5 < E 1 0 ♦ 3 , 3 X ) ) 


AU3072 

24 

180 

FORMAT <23X,9HRMS VALUE, 9 X,c (E l0.3*3X) ) 


AU3072 

25 

125 

FORMAT (/,23X,20HMEAN FILM THICKNESS ,5X, 


OC1872 

2 

1 

22HPR0B. OF SOLID CONTACT , / , 3 1 X , 6H < INCH) ,/ ) 


AU3072 

27 

130 

FORMAT (28X, FI 0.6, 16X, El 0.3) 


AU3072 

28 

195 

FORMAT <///5X, 10HMATR1X M , 36X , 1 OHM AT R I X C , 36X , 1 OHmATR IX K > 

VIBRAN 

94 

f90 

FORMAT (5x,3(E10.3,3X> ,7x,3(E10.3,3X> ,7X»3(E10.3,3X> ) 


VIBRAN 

95 


END 


VIBRAN 

96 
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SUBROUTINE READIT (TITlE»EN6m»EN6C»EN6k»NRAn6»SpEC»WInIT»WFIn»WDELT READIT 2 



1 ,FAC ) 

READIT 

3 


INTEGER SPEC 

READIT 

4 


DIMENSION ENGM(3),ENGC(3>,ENGK(3),WIMT(5>,WFIN<5>,WDELT<5* 

READIT 

5 


DIMENSION tItLE(G) 

READIT 

6 


COMMON/FT/ FLMTH(IO) »NEILMt 

AU3072 

29 


READ (5. 115) (TITLE (I) ,1=1.8) 

READIT 

7 


READ (5 , 100 > FAC 

READIT 

8 


DO 10 1=1,3 

READIT 

9 

10 

READ (5, 100 ) ENGM < I ) »ENGC ( I ) ,ENGK ( I ) 

READIT 

10 


READ (5,110) nrang.spec.nfilmt 

AU3072 

30 


DO 20 1=1 ,NRANG 

READIT 

12 

20 

REAO<5,100)WINIT(I) ,WFIN( I) ,WDELT ( I ) 

READIT 

13 


IF(nFIlmT.EO.O) go TO 500 

AU3072 

31 


READ (5, 100) (FLMTH ( I ) , 1=1 ,NFILMT) 

AU3072 

32 

100 

FORMAT (8E10. 3) 

READIT 

1A 

110 

FORMAT (515) 

READIT 

l5 

115 

FORMAT ( 8 A 1 0 ) 

READIT 

1 6 

500 

RETURN 

AU3072 

33 


END 

READIT 

18 
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SUBROUTINE WRIT1T ( T I TlE ,EnGm *EnGC , EnGk ,NRAnG , SPEC, 

WInIT.WFJn.WDElT 

WRITIT 

2 


1,FAC ) 


WRIT it 

3 


INTEGER SPEC 


WRITIT 

A 


DIMENSION ENGM (3) ,ENGC<3) ,ENGK (3) , WIN IT <5) ,WFIN(5) 

,wdelt(5* 

WRITIT 

5 


DIMENSION TITLE(B) 


wRItit 

6 


1=3 


WRITIT 

7 


WRITE (6*100) (TITLE (I) ,1 = 1 ,fi> 


WRITIT 

6 


WRITE(6,105)FAC 


WRITIT 

9 


WRITE<6»110)SPEC 


WRITIT 

10 


WRI TE (6, 120 ) 


WRITIT 

11 


WRITE (6, 130) 


WRITIT 

12 


DO 10 J=1 » I 


WRITIT 

13 

10 

WRITE(6,1A0) J,ENGM(J) ,ENGC(J) ,ENGk(J> 


WRITIT 

1A 


WRITE(6,150)nRAnG 


WRITIT 

15 


WRITE (6, 160 ) 


WRITIT 

16 


DO 20 1=1 , NRANG 


WRITIT 

17 

20 

WRITE (6,170) I.WINIT (I) ,WFIN(I) ,WDELT(I> 


WRITIT 

16 


return 


WRITIT 

19 

100 

FORMAT (1H1//30X, BA 10//) 


WRITIT 

20 

105 

FORMAT(A0X,27HMULTIPLIER for INPUT PSD = ,E10.3 /> 


WRITIT 

21 

1 1 o 

FORMAT (a0X,60HINPUT VIRRATION POWER SPECTRAL DENSITY SPECIFICATION 

WRITIT 

22 

1 CODE lS»I5»lH./> 


WRITIT 

23 

120 

FORMAT (AOX,65hROTOR BEARING SYSTEM PROPERTIES APE 

SPECIFItU by the 

WRITIT 

2A 


1 FOLLOWING /AOX, 15HSETS OF VALUES /) 


WRITIT 

25 

130 

FORMAT <ASX,50h MASS WEIGHT DAMPING 

SPRING / 

WRITIT 

26 


1 ASx,50HSPRlNG OF MASS (LBS*SEC 

stiffness / 

WRITIT 

27 


2 A5X , 50h NO (LBS) PER INCH) 

(LBS/INCh) ) 

WRITIT 

28 

1 AO 

FORMAT (ASX, 15,2 (5X,F 10.3) ,5X, El 0.3) 


WRITIT 

29 

150 

FORMAT(//AOX, 13HTHE FoLLoW I nG, 15, 30H FREQUENCY 

RANGES WILL be r 

WRITIT 

30 


1UN /) 


WRITIT 

31 

160 

FORMAT (ASX ,50HRANGE INITIAL FINAL 

frequency/ 

WRITIT 

32 


1 A5X , 50H NO FREQUENCY FREQUENCY 

INCREMENT) 

WRITIT 

33 

170 

FORMAT (A5X, 15,2 (5X ,Fl 0.2) »5X,F10.3) 


WRITIT 

3A 


END 


WRITIT 

35 
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SUBROUTINE SPEC1 <F, SCALE) 

SPEC1 

2 


SPEC1 

3 

SNAP8 SPEC 

SPEC1 

4 

NASA LEWIS SPEC 417-2. REV C, 1 JUNE69 

SPEC1 

5 

PARAGRAPH 3. 5. 1.2. 2 

SPEC1 

6 

EGS SYSTEM 50PERATING) 

SPEC1 

7 


SPEC1 

8 

IF (F .GT. 100. .AND. F .LT. 2000.) GO TO 5090 

SPEC1 

9 

IF (F .GT.2000. .OR. F.LT. 19.) GO TO 5110 

SPEC1 

10 

F X=F 

SPEC1 

11 

F2=l 00. 

SPEC1 

12 

PSD2=. 015 

SPEC1 

13 

R = 3. 

SPEC1 

14 

SCALE=pS02«(FX/F2 ) o*(r/3.) 

SpECl 

15 

RETU PN 

SPEC1 

16 

5090 SCALE=.015 

SPEC1 

17 

RETURN 

SPEC1 

18 

5110 SCALE=0- 

SPECl 

19 

RETURN 

SPEC1 

20 

END 

SPEC1 

21 
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12 

SUBROUTINE SPEC2 < F . SC ALE > 

SPEC2 

2 

SCALE= (2. *3. 14159*F) «*4/386.**2 

SPEC2 

3 

RETURN 

SPEC2 

4 

End. 

SPEC2 

5 
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13 

SUBROUTINE SPEC3(F. SCALE) 

SPEC3 

2 

scalf=.oi 

SPEC3 

3 

return 

SPEC3 

4 

END 

SPEC3 

5 
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SUBROUTINE MlNV ( A ,N »D *L »M > 

INVERT A MATRIX 

USAGE 

CALL MlNV <A,N,D,L»M> 

DESCRIPTION OF PARAMETERS „ Dw 

A - INPUT MATRIX, DISTROYED IN COMPUTATION AND REPLACED BY 

resultant inverse. 

N - ORDER OF MATRIX A 
D - RESULTANT DETERMINANT 
L - WORK VECTOR OF LENGTH N 
M - WORK VECTOR OF LENGTH N 

REMARKS 

MATRIX A MUST BE A GENERAL MATRIX 
METHOD 

THE standard gauss-jordan method is used, the determinant 

IS ALSO CALCULATED. A DETERMINANT OF ZERO INDICATES THAT 
THE MATRIX IS SINGULAR. 

DIMENSION A(3,1),L(1),M(1) 

COMPLEX A,BIGA,HOLO,D 
SEARCH FOR LARGEST ELEMENT 
D = 1.0 
NK=-N 

00 80 K=1,N 
NK=NK +N 
L ( K ) =K 
M(K)=K 
KK=NK+K 
EIGA=A(KK) 

DO 20 J=K,N 
IZ=N*( J-l) 

DO 20 I=K,N 
IJ=IZ*I 

10 IF(CABS(BIGA)-CABS(A(IJ) ) ) 15*20,20 
15 eiGA=A(IJ) 

L(K)=I 
M(K)=J 
20 CONTINUE 

INTERCHANGE ROWS 
J=L (K) 

IF(J-K) 35,35,25 
25 K I =K-N 

00 30 1=1, N 
K I=K I *N 
HOLD=-A(KI) 

JI=KI-K*J 
A (KI)=A(JI) 

30 A ( JI ) =HOLO 

INTERCHANGE COLUMNS 
35 I=M(K) 

IF(I-K) 45*45,38 
38 JP=N*(I-D 
DO 40 J=I,N 
JK=NK+J 
Ol=JP*J 
HOLD=-A(JK) 

A(JK)=A(JI) 


FEBI73 
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MlNV 

2 

MINV 

3 

MINV 

4 

MINV 

5 

MINV 

6 

MINV 

7 

MINV 

8 

MINV 

9 

MINV 

10 

MINV 

n 

MINV 

12 

MINV 

13 

MINV 

14 

MINV 

15 

MINV 

16 

MINV 

17 

MINV 

18 

MINV 

19 

MINV 

20 

MINV 

21 

M1 N V 

22 

MINV 

23 

MINV 

24 

MINV 

25 

MINV 

26 

MINV 

27 

MINV 

28 

MINV 

29 

MINV 

30 

MINV 

31 

MINV 

32 

MINV 

33 

MINV 

34 

MINV 

35 

MINV 

36 

MINV 

37 

MINV 

38 

MINV 

39 

MINV 

40 

MINV 

41 

MINV 

42 

MINV 

43 

MINV 

44 

MINV 

45 

MINV 

46 

MINV 

47 

MINV 

48 

MINV 

49 

MINV 

50 

MINV 

51 

MINV 

52 

MINV 

53 

MINV 

54 

MINV 

55 

MINV 

56 

MINV 

57 

MINV 

58 

MINV 

59 

MINV 

60 

/MlNV 

61 

MINV 

62 
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40 A ( JI ) -HOLD 

oivioe column by minus pivot (value of pivot element is 

CONTAINED IN BIGA) 

45 IP(CABS<BIGA) ) 48,46*48 

46 0=0.0 
RETURN 

48 DO 55 1=1, N 

IF(I-K) 50*55,50 
50 IK=NK+ I 

A(IK)=A(IK)/(-BlGA) 

55 CONTINUE 

REDUCE MATRIX 
DO 65 1 = 1, N 
I K=NK ♦ I 
HOLD=A ( IK) 

IJ=I-N 
DO 65 J=1,N 
I J= I J + N 

IE(I-K) 60,65,60 
60 IE(J-K) 62,65,62 
62 K J= I J- I +K 

A(IJ)=HOLD*A(KJ)+A(Ij> 

65 CONTINUE 

DIVIOE ROW BY PIVOT 
KJ=K-N 
DO 75 J=1*N 
K J=K J+N 

IE(J-K) 70,75,70 
70 A(KJ)=A(KJ)/BIGA 
75 CONTINUE 

PRODUCT OF PIVOJS 
D=D*BIGA 

REPLACE PIVOT BY RECIPROCAL 
A(KK)=1.0/BIGA 
80 CONTINUE 

FINAL ROW ANO COLUMN INTERCHANGE 
K=N 

100 K= ( K — 1 ) 

IF(K) 150,150,105 
105 I =L ( K ) 

IF(I-K) 120,120,108 
108 JQ=N*(K-1) 

JR=N*(I-1) 

DO HO U=1*N 
JK= JO + J 
HOLD=A< JK) 

JI=JR+J 

A(JK)=-A(JI) 

110 A ( JI ) =HOlD 
120 J=M(K) 

IF(J-K) 100,100,125 
125 KI=K-N 

DO 130 1=1, N 
K I =K I +N 
H0LD=A (KI ) 

JI=KI-K+J 
A (KI ) =-A < JI ) 

130 A ( JI ) =HOLD 
GO TO 100 
150 RETURN 
END 


M1NV 

63 

MINV 

64 

MINV 

65 

MINV 

66 

MINV 

67 

MINV 

68 

MINV 

69 

MINV 

70 

MINV 

71 

MINV 

72 

MINV 

73 

MlNV 

74 

MINV 

75 

MINV 

76 

MINV 

77 

MINV 

78 

MINV 

79 

MINV 

80 

MINV 

81 

MINv 

82 

MINv 

83 

MINv 

84 

MINV 

85 

MINV 

86 

MINV 

87 

MINV 

88 

MINV 

89 

MINV 

90 

MINV 

91 

MINV 

92 

MINV 

93 

MINV 

94 

MINV 

95 

MINV 

96 

MINV 

97 

MINV 

98 

MINV 

99 

MINV 

100 

MINV 

101 

MINV 

102 

MINV 

103 

MINV 

104 

MINV 

105 

MINV 

106 

MINV 

107 

MINV 

108 

MINV 

109 

MINV 

110 

MINV 

111 

MINV 

112 

MINV 

113 

MINV 

m 

MINV 

115 

MINV 

116 

MINV 

117 

MINV 

118 

MINV 

119 

MINV 

120 

MINV 

121 

MINV 

122 

MINV 

123 
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SUBROUTINE MaTRX2(M*C,K*ENGM,ENGC»£NGK) 

MATRX2 

2 

DIMENSION M (3,3) »C (3»3) *K (3*3) *D (3) »W (3) »Z£TA (3) *MU (3) 

MATRX2 

3 

DIMENSION ENGM ( 3 ) ,£NGC<3) »£NGK (3) 

MATRX2 

4 

REAL M*K ,MU 

MATRX2 

5 

DATA N/3/ 

MATRX2 

6 

00 1090 1 = 1, N 

MATRX2 

7 

ta(l)=SQRT (£NGK ( I ) / (ENGM ( I ) /386. ) > 

MATRX2 

8 

2ETA < I ) =.5« (ENGC < I ) / (SORT (ENGK < I ) *ENGM < I ) /386. ) ) ) 

MATRX2 

9 

1090 CONTINUE 

MATRX2 

10 

CO 1135 I = 1 * N 

MATRX2 

11 

CO 1135 J=1,N 

MATRX2 

12 

M ( I , J) =0 » 

MATRX2 

13 

IF(I.EQ.J) M(I,J)=1. 

MATRX2 

l4 

1135 CONTINUE 

MATRX2 

15 

MU ( 1 ) =1 • 

MATRX2 

16 

MU(2)=ENGM(2)/ENGM<1) 

MATRX2 

17 

MU ( 3 ) =ENGM <3) /ENGM (2) 

MATRX2 

10 

C(1»1)-ZETA(1)»«M1> 

MATRX2 

19 

C <1 *2) =-ZETa (2) *W (2)«MU (2) 

MATRX2 

20 

C (1,3)=0. 

MATRX2 

21 

C(2,1)=-ZeTA(1)«w(1) 

MATRX2 

22 

C(2,2)=(1.+MU<2) )*2ETA(2)*w(2) 

MATRX2 

23 

C <2.3)=-2ETA (3)*W(3)*Mu(3) 

MATRX2 

24 

C (3,1) =0.0 

MATRX2 

25 

C (3,2)=-ZETA(2)*W(2> 

MATRX2 

26 

C(3,3)=(1.*MU(3) )*ZETA(3)*to(3> 

MATRX2 

27 

K ( 1 ,1 ) = W (1 ) «*2 

MATRX2 

28 

K < 1 ,2) =-W (2) **2*MU (2) 

MATRX2 

29 

K ( 1 , 3) =0 • 

MATRX2 

30 

K (2,1) =-w ( 1 ) **2 

MATRX2 

31 

K (2,2)=W(2)«#2*(1 ,*MU(2) ) 

MATRX2 

32 

K (2,3) =~W (3 ) **2*MU ( 3 ) 

MATRX2 

33 

K ( 3, 1 ) =0 • 0 

MATRX2 

34 

M3»2)=-W(2)**2 

MATRX2 

35 

K (3,3) =w (3) **2*(1 .♦MU (3) ) 

MATRX2 

36 

RETURN 

MATRX2 

37 

END 

MATRX2 

38 
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16 

SUBROUTINE CAlC(M»C*K»W,PSD> 

CALC 

2 

DIMENSION M(3,3> »C(3»3> »K(3,3> *PsD(3) , 2 ( 3, 3) ,H (3, 3> ,DUMy1 <3> • 

CALC 

3 

1 DUMy2(3) 

CALC 

4 

COMPLEX EYE,2*DUMYl»DUMY2*H*DET2 

CALC 

5 

REAL M,K 

CALC 

6 

DATA N/3/ 

CALC 

7 

£ YE=CMPLX ( 0 • , 1 « ) 

CALC 

8 

DO 60 J= 1 , N 

CALC 

9 

DO 60 1=1, N 

CALC 

10 

60 2 (I,J) = (-to«»w*M(I,J)*K(I,j))+£YE*w*C(I,J) 

CALC 

11 

CALL MINV (2,N,DETZ,D0MY1 ,DUMY2) 

CALC 

12 

DO 110 1 = 1, N 

CALC 

13 

DO 110 J=1»N 

CALC 

14 

h<I,J)=2<I,J) 

CALC 

15 

PSD(i)=(CAeS{ H (l*l)))**2 

CALC 

16 

110 CONTINUE 

CAlC 

17 

RETURN 

CALC 

18 

END 

CALC 

19 




FEB173 

17 


SUBROUTINE TLU(A,B,C,D,N> 

AU3072 

34 


DIMENSION C < 1 > ,0 (1) 

AU3072 

35 


IF (N-l ) 1 ,2 ,3 

AU3072 

36 

1 

B=0 * 

AU3072 

37 


C-0 TO 100 

AU3072 

38 

2 

e=D(l) 

AU3072 

39 


C-0 TO 100 

AU3072 

40 

3 

ML = 1 

AU3072 

41 


MU=N 

AU3072 

42 

8 

IF (MU-ML-1 ) 15,15,9 

AU3072 

43 

9 

M= (MU + ML) /2 

AU3072 

44 


IF(C(1)-C(2) ) 11,2,10 

AU3072 

45 

10 

1F(C(M)-A) 13,12, 14 

AU3072 

46 

11 

IF(A-C(M) ) 13,12,14 

AU3072 

47 

12 

B=D ( M ) 

AU3072 

48 


C-0 TO 100 

AU3072 

49 

13 

MU=M 

AU3072 

50 


C-0 TO 8 

AU3072 

51 

14 

ML=M 

AU3072 

52 


GO TO 8 

AU3072 

53 

15 

R=D ( ML ) ♦ (D (Mu)-D (ML) >*( ( A-C(ML) ) / ( C (MU) -C ( ML ) ) ) 

AU3072 

54 

100 

RETURN 

AU3072 

55 


END 

AU3072 

56 
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SAMPLE TEST HARE IHPOT 

1 

axial random vibration response - nominal case 10-10-72 


0.01 

96. 

0.02 

300000. 

7. 

0.02 

250000. 

21. 

84. 

113500. 

3 

1 5 


20. 

200. 

2. 

200. 

400. 

5. 

400. 

2000. 

20. 


.0004 .0005 .0006 .0007 .0008 


SAMPLE TEST CASE OUTPtlT 


AXIAL RANDOM VIBRATION RESPONSE - NOMINAL CASE 10-10-72 


MULTIPLIER FOR INPUT PSD = 1.000E-02 

INPUT VIBRATION POWER SPECTRAL DENSITY SPECIFICATION CODE IS 1. 

ROTOR BEARING SYSTEM PROPERTIES ARE SPECIFIED BY THE FOLLOWING 
SETS OF VALUES 


MASS WEIGHT 

SPRING OF MASS 

NO (LBS) 

1 96.000 

2 7.000 

3 21.000 


DAMPING 
(LBS'SEC 
PER INCH) 
.020 
.020 
BA. 000 


SPRING 
STIFFNESS 
(LBS/INCH) 
3.000E»05 
2.500E*05 
1 • 135E-05 


THE FOLLOWING 

RANGE 

NO 

1 

2 

3 


3 FREQUENCY 

INITIAL 

FREQUENCY 

20.00 

200.00 

AOO.OO 


RANGES WILL 

FINAL 

FREQUENCY 

200.00 

400.00 

2000.00 


BE RUN 

FREQUENCY 

Increment 

2.000 

5.000 

20.000 


MATRIX M 
1 • 000E+00 

0 . 

0 . 

MATRIX C 
4.021E-02 

-4.021E-02 

0 . 

MATRIX K 
1.206E*06 

-1.005E*06 

0 . 

0 . 

1 d'OOOE + OO 

0 . 

-4*021E-02 

5.916E-01 

-2.316E*03 

-I.206E*06 

1.479E*07 

-6.259E*06 

0 . 

0 . 

1.000E*00 

0 . 

-5.514E-01 

3.088E»03 

0. 

-I.379E»07 

8.345E*04 



FREQUENCY RANGE NUMBER 1 

COVERS THE RANGE FROM 20.00 TO 200.00 
IN INCREMENTS OF 2.00 


FREQUENCY 

PSO OF Yl 

PSD OF Y2 

PSD OF Y3 

PSD INPUT 

PSO INPUT 

(HZ) 

UN* IN/HZ) 

( IN® IN/HZ i 

(IN*IN/HZ) 

<IN«IN/HZ> 

1G*G/HZ> 

20.000 

S.328E-12 

3.9756-13 

1.087E-12 

1.792E-08 

3.000E-05 

22.000 

5.910E-12 

4. 4246-13 

1.210E-12 

V.347E-08 

. 3.300E-0S 

24.000 

6.506E- 1 2 

4.888E-13 

1.337E-12 

1.037E-08 

3.600E-05 

26.000 

7.119E-12 

S.370E-13 

1.470E-12 

8.159E-09 

3.900E-0S 

28.000 

7.750E-12 

S.871E-13 

1.607E-12 

6.S32E-09 

4.200E-05 

30.000 

8 .40 IE- 12 

6.394E-13 

1.751E-12 

5.311E-09 

4.SOOE-OS 

32.000 

9.074E-12 

6.9426-13 

1.902E-12 

4.376E-09 

4.800E-05 

34.000 

9.772E-12 

7.5166-13 

2.061E-12 

3.64BE-09 

S.100E-05 

36.000 

1 . 050E-1 1 

8.119E-13 

2.227E-12 

3.074E-09 

5.400E-05 

38.000 

1.125E-11 

8.7S5E-13 

2.403E-12 

2.613E-09 

5.700E-05 

40.000 

1.204E-11 

9.4276-13 

2.589E-12 

2.241E-09 

6.000E-05 

42.000 

1.286E-U 

1.014E-12 

2.787E-12 

1.9366-09 

6.300E-05 

44.000 

1 .372E-11 

1.090E-12 

2.997E-12 

1.683E-09 

6.6006-05 

46.000 

1.463E-11 

1.1706-12 

3.220E-12 

1.473E-09 

6.9006-05 

48.000 

1.550E-11 

1.256E-12 

3.4S9E-12 

1.297E-09 

7.20OE-OS 

50.000 

1.658E-11 

1.348E-12 

3.715E-12 

1.147E-09 

7.500E-05 

' 52.000 

1.764E-11 

1.446E-12 

3.990E-12 

1.020E-09 

7.800E-05 

54.000 

1.876E-11 

1.5526-12 

4.285E-1 2 

9.107E-10 

8.1006-05 

56.000 

1.995E-11 

1.666E-12 

4.603E-12 

8.1666-10 

8.400E-05 

58.000 

2.122E-11 

1.789E-12 

4.947E-12 

7.350E-10 

8.700E-05 

60.000 

2.257E-1 1 

1.921E-I2 

5.320E-12 

6.6396-10 

9.0006*05 

62.000 

2.401E-11 

2.066E-12 

5.725E-12 

6.017E-10 

9.300E-05 

64.000 

2.556E-1 1 

2.2236-12 

6.166E-12 

5.470E-10 

9.600E-05 

66.000 

2.722E-11 

2.3946-12 

6.647E-12 

4.9886-10 

9.900E-05 

68.000 

2.900E-11 

2.581E-12 

7 . 174E-12 

4.56IE-10 

1.020E-04 

70.000 

3.093E-11 

2.7866-12 

7.752E-12 

4.1812*10 

1.050E-04 

72.000 

3.302E-1 1 

3.011E-12 

8.388E-I2 

3.8426-10 

1.080E-04 

74.000 

3.529E-11 

3.2S9E-12 

9.090E-12 

3.539E-10 

1.1 10E-04 

76.000 

3.776E-11 

3.5346-12 

9.868E-12 

3.267E-1 0 

1.140E-04 

78.000 

4.04SE-11 

3.8396-12 

1.073E-1I 

3.022E-I0 

1.170E-04 

80.000 

4,341 E— 1 1 

4.178E-12 

1.170E-11 

2.B01E-10 

1.200E-04 

82.000 

4.666E-1 1 

4.557E-12 

1.277E-11 

2.601E-10 

1.230E-04 

04.000 

5. 025E-1 1 

4.983E-12 

1.398E-11 

2.419E-10 

1.260E-04 

86.000 

5.423E-11 

5.461E-12 

1.53SE-11 

2.255E-10 

1.290E-04 

88.000 

5.065E-11 

6.0036-12 

1.689E-11 

2.104E-10 

1.320E-04 

90.000 

6.360E-1 1 

6.6186-12 

1.065E-11 

1.967E-10 

1.350E-04 

92.000 

6.916E-1 1 

7.3206-12 

2.065E-11 

1.842E-10 

1 • 380E-04 

94.000 

7.544E-11 

8.1266-12 

2.296E-11 

1.726E-10 

1.410E-04 

96.000 

8.257E-II 

9.0566-12 

2.562E-II 

1.62JE-IO 

I .4406-04 

98.000 

9. 072E-1 1 

1.014E-li 

2.072E-11 

1.524E-10 

1.4706-04 

100.000 

1.001E-10 

1.140E-li 

3.235E-11 

1.434E-10 

1. 5006-04 

102.000 

1.088E-10 

1.263E-li 

3.591E-11 

1.325E-10 

1.5006-04 

104.000 

1 .189E-10 

1.409E-11 

4.010E-11 

1.2266-10 

1.5006-04 

106.000 

I.307E-10 

1.5826-11 

4.510E-11 

1.1366-10 

1. 5006-04 

108.000 

1.447Et10 

1.7906-11 

5.110E-1I 

1.054E-10 

' 1.5006-04 

110.000 

1.615E-10 

2.042E-I1 

5.840E-11 

9.794E-11 

1.5006-04 

112.000 

1.8I9E-10 

2.352E -ll 

6.737E-H 

9.H3E-11 

1.5006-04 

114.000 

2.069E-10 

2.738E -li 

7.857E-I1 

8.490E-11 

1.5006-04 

116.000 

2.382E-10 

3.2286-U 

9.278E-11 . 

7.920E-11 

1.5006-04 

118.000 

2.780E-10 

3.860E-11 

1*111 E— 1 0 

7.396E-11 

1 .5006-04 

120.000 

3.297E- 1 0 

4.6956-11 

1.354E-10 

&.9I5E-11 

1.5006-04 

122.000 

3.986E-10 

5.826E-11 

1.683E-10 

6.473E-11 

1.500E-04 

124.000 

4.935E-10 

7.407E-U 

2.144E-10 

6.0656-11 

1.5006-04 

126.000 

6.289E-10 

9.7016-11 

2.813E-10 

S.689E-11 

1.5006-04 

128.000 

8.313E-10 

1.319E-10 

3.830E-10 

5.342E-11 

1.5006-04 

130.000 

1.151E-09 

1.879E-10 

5.469E-10 

S.021E-11 

1.5006-04 

132.000 

1.692E-09 

2.845E-10 

8.293E-10 

4.723E-11 

1.5006-04 

134.000 

2.674E-09 

4.6346-10 

1.353E-09 

4.448E-11 

1.5006-04 

136.000 

4.513E-09 

8.088E-10 

2.361E-09 

4,1 92E-1 I 

1.5006-04 

130.000 

7.241E-09 

1.336E-09 

3.918E-09 

3.954E-11 

1.5006-04 

140.000 

7.807E-09 

1.489E-09 

4.373E-09 

3.733E-11 

1.S00E-04 

142.000 

5. 096E-09 

1.005E-09 

2.957E-09 

3.5276-11 

1.500E-04 

144.000 

2.877E-09 

5.8736-10 

1.732E-09 

3.335E-11 

1.5006-04 

146.000 

1.703E-09 

3.601E-10 

1.064E-09 

3.156E-11 

1.S006-04 

148.000 

1 •O06E-O9 

2.382E-10 

7.050E-10 

2.989E-11 

1.500E-04 

150.000 

7 .39 IE- 10 

1.681E-10 

4.986E-I0 

2.833E-11 

1.5006-04 

152.000 

5.289E-1 0 

1.249E-10 

3.713E-10 

2.6866-11 

1.500E-04 

154.000 

3.937E-10 

9.669E-U 

2.879E-10 

2.550E-11 

1.5006-04 

156.000 

3.024E-10 

7.727E-11 

2.306E-10 

2.42)6-11 

1.5006-04 

150.000 

2.383E-10 

6.339E-H 

1.096E-1O 

2.301E-11 

1.5006-04 

160.000 

1.916E-10 

5.314E-11 

1.592E-10 

2. 1886-1 1 

1.500E-04 

162.000 

1.568E-10 

4.537E-U 

1.362E-10 

2.082E-11 

1.500E-04 

164.000 

1.302E-10 

3.933E-1 1 

1.183E-10 

1.982E-11 

1.500E-04 

166.000 

1.095E-10 

3. 4566-11 

1.042E-10 

1.888E-U 

l.SOOE-04 

160.000 

9.311 E- II - 

3.073E-11 

9.282E-22 

l.BOOE-ll 

1.5006-04 

170.000 

7.994E-11 

2.7606-U 

8.355E-11 

1.717E-11 

1.500E-04 

172.000 

6.922E-1 1 

2.502E-11 

7.591E-11 

1.638E-11 

1.5006-04 

174.000 

6.041E-11 

2.288E-11 

6.954E-11 

1.5646-11 

1.5006-04 

176.000 

5.309E-11 

2. 107E-1 1 

6.419E-11 

1.49SE-11 

1.500E-04 

i7e.ooo 

4.698E-11 

1. 9546-11 

5.966E-11 

1.428E-11 

1.5006-04 

180.000 

4. 182E-1 1 

1.824E-11 

5.580E-11 

1.3666-11 

1.500E-04 

182.000 

3.746E-1 1 

1.7136-11 

5.250E-11 

1.307E-11 

1. 5006-04 

184.000 

3.374E-11 

1.6176-U 

4.966E-11 

1.2S1E-11 

1.S00E-04 

186.000 

3.056E-11 

I. 5346- U 

4.721E-I1 

I.198E-11 

1.5006-04 

188.000 

2.784E-11 

1.4626-U 

4.510E-11 

1.148E-11 

1.5006-04 

190.000 

2.551E-U 

1.4006-11 

4.327E-11 

1.100E-11 

1.5006-04 

192.000 

2.350E-1 1 

1.3456-11 

4.168E-11 

1.055E-11 

1.500E-04 

194.000 

2.178E-11 

1.298E-11 

4.030E-11 

1.012E-11 

1.S00E-04 

196.000 

2.031E-1 1 

1.257E-I1 

3.91 IE-11 

9.717E-12 

1.500E-04 

198.000 

1.906E-11 

1.221E-11 

3.808E-11 

9.330E-1Z 

1.500E-04 

200.000 

1.801E-11 

1.190E-U 

3.719E-11 

8.962E-12 

1.500E-04 

INTEGRATED value 

8.972E-08 

1. 6996-08 

4 . 993E-08 

I.772E-07 

2.2206-02 

total int. value 

8.972E-08 

1.6996-08 

4.993E-08 

I.772E-07 

2.2206-02 

RMS VALUE 

2.995E-04 

1.3046-04 

2.235E-04 

4.209E-04 

1.4906-01 
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FREQUENCY RANGE 

NUMBER 2 





COVERS THE RANGE 

FROM 200. 

00 TO 400 

.00 



IN INCREMENTS OF 

5.00 





FREQUENCY 

PSD OF Y1 

PSD OF Y2 

PSD Of Y3 

PSD INPUT 

PSD INPUT 

(HZ) 

<IN*IN/HZ) 

( IN*IN/HZ ) 

(IN*IN/HZ) 

(INMN/HZ) 

(G*G/HZ) 

200.000 

1.801E-11 

1 . 190E-1 1 

3.719E-11 

8.962E-12 

1.500E-04 

205.000 

1.611E-11 

1.128E-U 

3.543E-11 

8.120E-12 

1.500E-04 

210.000 

1.507E-11 

1.082E-11 

3.416E-11 

7.373E-12 

1.500E-04 

215.000 

1.473E-11 

1.043E-H 

3.309E-11 

6.71 IE-12 

1.500E-04 

220.000 

1.493E-11 

1.002E-U 

3.194E-11 

6.12iE-12 

1.500E-04 

225.000 

1.547E-11 

9.493E-12 

3.041E-11 

5.595E-12 

1.500E-04 

230.000 

1.608E-11 

8.768E-12 

2.821E-11 

5.124E-12 

1.500E-04 

235.000 

1.646E-11 

7 .816E-12 

2.527E-11 

4.702E-12 

1.500E-04 

240.000 

1.639E-11 

6.689E-12 

2.172E-11 

4.322E-12 

1.500E-04 

245.000 

1.579E-11 

5.501E-12 

1.793E-11 

3.980E-12 

1.500E-04 

250.000 

1.473E-11 

4.375E-12 

1.432E-11 

3.671E-12 

1.500E-04 

255.000 

1.340E-11 

3.398E-12 

1.116E-11 

3.391E-12 

1.500E-04 

260.000 

1.199E-11 

2.604E-12 

8.582E-12 

3.138E-12 

1.500E-04 

265.000 

1.062E-11 

1.985E-12 

6.562E-12 

2.908E-12 

1.500E-04 

270.000 

9.357E-12 

1.514E-12 

5.020E-12 

2.698E-12 

1.500E-04 

275.000 

8.240E-12 

1.161E-12 

3.857E-12 

2.507E-12 

1.500E-04 

280.000 

7.265E-12 

8.960E-13 

2.983E-12 

2.333E-12 

1.500E-04 

285.000 

6.423E-12 

6.975E-13 

2.327E-12 

2.174E-12 

1.500E-04 

290.000 

5.698E-12 

5.480E-13 

1.830E-12 

2.027E-12 

1.500E-04 

295.000 

5.075E-12 

4 . 345E-1 3 

1.452E-12 

1.893E-12 

1.500E-04 

300.000 

4.537E-12 

3.477E-1 3 

1 • 163E-12 

1.770E-12 

1.500E-04 

305.000 

4.072E-12 

2.807E-13 

9.384E-13 

1.657E-12 

1.500E-04 

310.000 

3.668E-12 

2.286E-13 

7.635E-13 

1.553E-12 

1.500E-04 

315.000 

3.317E-12 

1.876E-13 

6.260E-13 

1.456E-12 

1.500E-04 

320.000 

3 • 009E-12 

1.552E-13 

5.169E-13 

1.368E-12 

1.500E-04 

325.000 

2.738E-12 

1.293E-13 

4.297E-13 

1.285E-12 

1.500E-04 

330.000 

2.499E-12 

1.085E-13 

3.595E-13 

1.209E-12 

1.500E-04 

335.000 

2.288E-12 

9.164E-14 

3.025E-13 

1.139E-12 

1.500E-04 

340.000 

2.099E-12 

7.789E-14 

2.560E-13 

1.073E-12 

1.500E-04 

345.000 

1.931E-12 

6.659E-14 

2.178E-13 

1.012E-12 

1.500E-04 

350.000 

1.781E-12 

5.726E-14 

1.863E-13 

9.556E-13 

1.500E-04 

355.000 

1.646E-12 

4.951E-14 

1.600E-13 

9.029E-13 

1.500E-04 

360.000 

1.524E-12 

4.302E-14 

1.381E-13 

8.538E-13 

1.500E-04 

365.000 

1.414E-12 

3.757E-14 

1.197E-13 

8.079E-13 

1.500E-04 

370.000 

1.314E-12 

3.297E-14 

1.041E-13 

7.651E-13 

1.500E-04 

375.000 

1.224E-12 

2.906E-14 

9.093E-14 

7.251E-13 

1.500E-04 

380.000 

1.141E-12 

2.573E-14 

7.970E-14 

6.877E-13 

1.500E-04 

385.000 

1.066E-12 

2.288E-14 

7.010E-14 

6.527E-13 

1.500E-04 

390.000 

9.972E-13 

2.042E-14 

6.186E-14 

6.199E-13 

1.500E-04 

395.000 

9.341E-13 

1.830E-14 

5 .476E”14 

5.891E-13 

1.500E-04 

400.000 

8.762E-13 

1.646E-14 

4.862E-14 

5.602E-13 

1.500E-04 

INTEGRATED VALUE 

1.432E-09 

5.320E-10 

1.708E-09 

5.232E-10 

3.000E-02 

TOTAL INT. VALUE 

9.115E-08 

1.752E-08 

5.164E-08 

1.777E-07 

5.220E-02 

RMS VALUE 

3.019E-04 

1.324E-04 

2.272E-04 

4.21SE-04 

2.285E-01 



frequency range number 3 

COVERS THE RANGE FROM 400.00 TO 2000.00 

IN INCREMENTS OF 20.00 


FREQUENCY 

PSD OF Y1 

PSD OF Y2 

PSD OF Y3 

PSD INPUT 

PSD INPUT 

(HZ) 

(IN*IN/HZ> 

(INMN/HZ1 

( IN* IN/HZ) 

(INMN/HZ) 

(6*6/H21 

400.000 

8.762E-13 

1.646E-1* 

4.862E-14 

S.602E-13 

l.SOOE-04 

420.000 

6.866E-13 

1.1166-14 

3.104E-14 

4.608E-13 

1.500E-04 

440.000 

5.474E-13 

7.9S4E-15 

2.060E-14 

3.826E-13 

1.500E-04 

460.000 

4.428E-13 

S.920E-1S 

1.412E-14 

3.203E-13 

l.SOOE-04 

480.000 

3.627E-13 

4.573E-18 

9.956E-1S 

2.701E-13 

l.SOOE-04 

500.000 

3.004E-J3 

3.646E-1S 

7.185E-15 

2.294E-J3 

1 .5006-04 

520.000 

2.51 3E-1 3 

2.984E-1S 

5.289E-15 

1.96XE-13 

l.SOOE-04 

540.000 

2.120E-13 

2.4936-15 

3.955E-15 

1.6B6E-13 

l.SOOE-04 

560.000 

1.803E-13 

2.1176-15 

2.995E-1S 

1.45BE-13 

l.SOOE-04 

580.000 

1.544E-13 

1.8176-15 

2.287E-I5 

1.267E-13 

l.SOOE-04 

600.000 

1.331E-13 

1.570E-15 

1.755E-IS 

1.1 066-13 

1.S00E-04 

620.000 

1.154E-13 

1.3586-15 

1.349E-1S 

9.705E-14 

l.SOOE-04 

640.000 

1 .005E-13 

1.172E-15 

1.034E-15 

8.547E-14 

l.SOOE-04 

660.000 

8.806E-14 

1.0066-15 

7.887E-16 

7.5S7E-14 

1.S00E-04 

680.000 

7.747E-14 

8.S54E-16 

5.971E-16 

6.7076-14 

1.500E-04 

700.000 

6.844E-14 

7.2016-16 

4.480E-16 

5.9726-14 

l.SOOE-04 

720.000 

6.069E-14 

S. 9976-16 

3.332E-16 

5.3366-14 

l.SOOE-04 

740.000 

5.402E-14 

4.943E-16 

2.457E-16 

4.782E-14 

l.SOOE-04 

760.000 

4.824E-14 

4.038E-16 

1.800E-16 

4.298E-14 

1.500E-04 

780.000 

4.322E-14 

3.276E-16 

1.311E-16 

3.874E-14 

l.SOOE-04 

800.000 

1.884E-14 

2.6446-16 

9.529E-17 

3.501E-14 

1.500E-04 

820.000 

3.500E-14 

2. 127E- 16 

6.918E-17 

3.172E-14 

l.SOOE-04 

840.000 

3.163E-14 

1.7096-16 

5.026E-17 

2.880E-14 

1.5006-04 

860.000 

2.865E-14 

1 .374E-16 

3.661E-17 

2.622E-14 

1. 5006-04 

880.000 

2.602E-14 

1.106E-16 

2.676E-17 

2.391E-14 

l.SOOE-04 

900.000 

2.369E-14 

8.927E-17 

1.965E-17 

2.186E-14 

1.500E-04 

920.000 

2.162E-14 

7.2276-17 

1.450E-I7 

2.002E-14 

1.5006-04 

940.000 

1.977E-14 

5.8726-1? 

1.076E-17 

1.837E-14 

l.SOOE-04 

960.000 

1.811E-14 

4.790E-17 

8.035E-18 

1.688E-14 

1.500E-04 

980.000 

1.663E-14 

3.9246-1? 

6.034E-18 

1.S5SE-14 

1.5006-04 

1000.000 

1.529E-14 

3.229E-1? 

4.559E-18 

1.434E-14 

1.500E-04 

1020.000 

1 .409E-14 

2.668E-1? 

3.466E-18 

1.325E-14 

1.500E-04 

1040.000 

1.301E-14 

2.214E-1? 

2.650E-18 

1.2266-14 

1.S00E-04 

1060.000 

1.202E-14 

1.845E-1? 

2.038E-18 

1.1 36E-14 

1.5006-04 

1080.000 

1.113E-14 

1.S43E-17 

1.576E-18 

1.0S4E-14 

l.SOOE-04 

1100.000 

1.032E-14 

1.2976-1? 

1 .226E-18 

9.794E-1S 

1.S00E-04 

1120.000 

9.587E-15 

1 . 094E-1 7 

9.587E-19 

9.1136-15 

1.500E-04 

1140.000 

8.916E-15 

9.2S7E-18 

7.535E-19 

8.490E-15 

1. 5006-04 

1160.000 

8.302E-15 

7.8656-18 

5.951E-19 

7.9206-15 

l.SOOE-04 

1180.000 

7.740E-15 

6.7066-18 

4.723E-19 

7.396E-1S 

l.SOOE-04 

1200.000 

7.226E-15 

S.736E-18 

3.766E-19 

6.915E-1S 

1.5006-04 

1220.000 

6.754E-15 

4.9236-18 

3.016E-19 

6.473E-15 

1.500E-04 

1240.000 

6.319E-15 

4.2386-18 

2.426E-19 

6.06SE-15 

1.S00E-04 

1260.000 

5.920E-15 

3.6606-18 

1.960E-19 

S.689E-15 

1.500E-04 

1280.000 

S.551E-15 

3.1696-18 

1.589E-19 

5.3426-15 

1.500E-04 

1300.000 

5.211E-15 

2.7526-18 

1.294E-19 

5.021E-15 

1.500E-04 

1320.000 

4.897E-15 

2.396E-18 

1.057E-19 

4.723E-15 

l.SOOE-04 

1340.000 

4.606E-15 

2.092E-18 

8.667E-20 

4.448E-15 

1 .500E-04 

1360.000 

4.336E-15 

1.8316-18 

7.131E-20 

4. 192E-15 

1.500E-04 

1380.000 

4.086E-15 

1.606E-18 

5.888E-20 

3.954E-15 

1.500E-04 

1400.000 

3.R54E-15 

1.412E-18 

4.876E-20 

3.733E-15 

l.SOOE-04 

1420.000 

3, 63BE-1 5 

1 .2442-18 

4.052E-20 

3.527E-15 

1.500E-04 

1440.000 

3.437E-15 

1.099E-18 

3.376E-20 

3.3356-15 

1.500E-04 

1460.000 

3.250E-15 

9.724E-19 

2.822E-20 

3. 156E-15 

1.500E-04 

1480.000 

3.075E-15 

8.623E-19 

2. 365E-20 

2.989E-15 

1.500E-04 

1500.000 

2.912E-15 

7.660E-19 

1.988E-20 

2.B33E-15 

1.500E-04 

1520.000 

2.760E-15 

6.8186-1$ 

1.675E-20 

2.686E-1S 

1.500E-04 

1540.000 

2.617E-15 

6.080E-19 

1.415E-20 

2.550E-15 

1.S00E-04 

1560.000 

2.484E-15 

5.4316-19 

1.198E-20 

2.421E-15 

1.500E-04 

1580.000 

2.359E-15 

4.8S9E-19 

1.017E-20 

2.301E-15 

l.SOOE-04 

1600.000 

2.242E-15 

4. 35SE- 1 9 

8.657E-21 

2.188E-15 

l.SOOE-04 

1620.000 

2.132E-15 

3.909E-19 

7.383E-21 

2.082E-1S 

1.500E-04 

1640.000 

2.029E-15 

3.S15E-19 

6.31 IE-21 

1.982E-15 

1.500E-04 

1660.000 

1.932E-15 

3.1656-19 

5.406E-21 

1.888E-15 

1.S00E-04 

1680.000 

1.840E-1S 

2 ■ 854E- 1 9 

4.640E-21 

1.800E-15 

1.S00E-04 

1700.000 

1 .754E-15 

2.577E-19 

3.991E-21 

1.717E-15 

l.SOOE-04 

1720.000 

1 .673E-15 

2 • 330E- 1 9 

3.439E-21 

1.638E-1S 

1.S00E-04 

1740.000 

1 .597E-1S 

2. UOE-19 

2.970E-21 

1 . 564E-15 

1.500E-04 

1760.000 

1 .525E-15 

1.913E-19 

2.569E-21 

1.495E-15 

1.500E-04 

1780.000 

1.457E-15 

1 . 736E-19 

2.226E-21 

1.428E-15 

1.500E-04 

1800.000 

1 .392E-15 

1.578E-19 

1 .933E-21 

1.366E-15 

1.500E-04 

1820.000 

1.332E-15 

1.436E-19 

1.681E-21 

1 .307E-15 

1.500E-04 

1840.000 

1.274E-15 

1.308E-I9 

1.464E-2I 

1.251E-1S 

1.S00E-04 

1860.000 

1 .220E-15 

1.1936-1$ 

1.278E-21 

1.198E-15 

1.500E-04 

1880.000 

1 .168E-15 

1 .089E-19 

1.117E-21 

1.148E-1S 

1.500E-04 

1900.000 

1.119E-15 

9.9S9E-20 

9.775E-22 

1 . 1 00E-1 5 

1.500E-04 

1920.000 

1 .073E-15 

9.113E-20 

8.570E-22 

1.05SE-15 

1.500E-04 

1940.000 

1 .029E-15 

8.347E-20 

7.524E-22 

1.0126-15 

1.500E-04 

1960.000 

9.874E-16 

7.6536-20 

6.616E-22 

9.71 7E-I6 

1 .5006-04 

1980.000 

9.470E-16 

7.024E-20 

5.825E-22 

9.330E-16 

1.S00E-04 

2000.000 

9.102E-16 

6.453E-20 

5.136E-22 

8. 9626-16 

1 .5006-04 

INTEGRATED VALUE 

9.702E-11 

1.2166-12 

2.579E-12 

7.428E-1 1 

2.400E-01 

total int. value 

9.1 25E-08 

1 . 753E-08 

5.164E-08 

1 .7786-07 

2.9226-01 

RMS VALUE 

3.021E-04 

1.324E-04 

2.272E-04 

4.2166-04 

S. 4066-01 

MEAN film THICKNESS 

PROS. 

OF SOL 10 CONTACT 




(INCH) 


.000400 3.936E-02 
.000500 I .3896-02 
.000600 4.180E-03 
.000700 1.059E-03 
.000800 2.257E-04 
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Accelerometer on , 
Itop of bru simulator 


INSTRUMETi^mON^ 

turret! 


VERTICA L SUPPORT! 
fixtureSHB 


ACCELEROMETER ON| 

TOP OF VERTICAL SUPPO RT 

fixtureHB 


VIBRATION TABLE 


Fig. 3 BRU Simulator In Vertical Support Fixture For Vibration Testing 


MTI-11039 











Fig. 4 BRU Simulator In Support Fixture For Transverse 
Vibration Testing 


MII-11265 


COMPRESSOR END UP 



21 

TURBINE END DOWN 


Fig. 5 Proximity Probe Locations on BRU Simulator Rotor 
and Bearings 


MTI-12289 



RANDOM VIBRATION INPUT 
LEVEL: 0.54 G RMS 

ONE MAJOR VERTICAL 
DIVISION EQUALS : 0.0002 INCH 

ONE MAJOR HORIZONTAL 
DIVISION EQUALS: 10 m SEC 



RANDOM VIBRATION INPUT 
LEVEL : 1.52 G RMS 

ONE MAJOR VERTICAL 
DIVISION EQUALS •• 0.0002 INCH 

ONE MAJOR HORIZONTAL 
DIVISION EQUALS » 10 m SEC 




AMPLITUDE 

PROBABILITY 


AMBIENT GAS= N 2 
AMBIENT PRESSURE = 35 PSIA 
ROTOR SPEED: 36,000 RPM 
HYDROSTATIC SUPPLY PRESSURE: 

THRUST BRG = — 

TURBINE JOURNAL BRG = — 
COMPRESSOR JOURNAL BEARING: 55 PSIA 
DISPLACEMENT PROBE NO. 13 


Fig. 6 Casing-To-Pad Leading Edge Amplitude-Time History For Flex-Mounted 
Turbine Journal Bearing Pad Under Externally-Imposed Shaped Random 
Vibrations According To NASA Spec 417-2-C-3.5 


MTI-14034 






UPPER TRACE: 


FLEX. MOUNTED PAD (D) 
TURPINF JOURNAL 

LOWER TRACE: FLEX. MOUNTED PAD (A) 

COMPRESSOR JOURNAL 

VERTICAL SCALE: 0.00025 IN./DIV. 

HORIZONTAL SCALE: 20 MSEC/DIV. 


UPPER TRACE: SOLID MOUNTED PAD (F) 

TURBINE JOURNAL 

LOWER TRACE: SOLID MOUNTED PAD (C) 

COMPRESSOR JOURNAL 

VERTICAL SCALE: 0.00025 IN./DIV. 

HORIZONTAL SCALE:. 20 MSEC/DIV. 


UPPER TRACE: COMPRESSOR JOURNAL 

FLEXURE PROBE (NO. 22) 

VERTICAL SCALE: 0.0005 IN./DIV. 

LOWER TRACE: THRUST BEARING 

GIMBAL PROBE (NO. 23) 

VERTICAL SCALE: 0.00075 IN./DIV. 
HORIZONTAL SCALE: 20 MSEC/DIV. 





RANDOM VIBRATION INPUT LEVEL: 0.95 G RMS 

AMBIENT CONDITIONS: N @ 14.7 PSIA 

2 

HYDROSTATIC SUPPLY PRESSURE: 150 PSIA 

ROTOR SPEED: 2500 RPM 



Fig. 7 


Bearing Amplitudes Under Shaped Random Vibrations According 
To NASA Spec 417-2-C-3.5 







ONE MAJOR DIVISION 
EQUALS 0.0005 IN. 


CAMERA SHUTTER SPEED: 1/50 SEC. 



RANDOM VIBRATION INPUT LEVEL: 1.2 G RMS 

AMBIENT CONDITIONS: N @ 14.7 PSIA 
2 

HYDROSTATIC SUPPLY PRESSURE: 150 PSIA 
ROTOR SPEED: 2500 RPM 

Fig. 8 Rotor Amplitudes (Casing-To-Shaf t) Under Shaped Random Vibrations 
According to NASA Spec 417-2-C-3.5 


MT-11733 










RANDOM VIBRATION INPUT 
LEVEL ■ 0.54 G RMS 

ONE MAJOR VERTICAL 
DIVISION EQUALS = 0.00025 INCH 

ONE MAJOR HORIZONTAL 
DIVISION EQUALS =' 10 msec 
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RANDOM VIBRATION INPUT 
LEVEL '- 1.52 G R M s 

ONE MAJOR VERTICAL 
DIVISION EQUALS : 0.00025 INCH 

ONE MAJOR HORIZONTAL 
DIVISION EQUALS = 10 msec 
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AMBIENT GAS = N 2 
AMBIENT PRESSURE : 35 PSIA 
ROTOR SPEED = 36,000 RPM 
HYDROSTATIC SUPPLY PRESSURE = 

THRUST BRG = — 

TURBINE JOURNAL BRG : — 
COMPRESSOR JOURNAL BEARING 55 PSIA 
DISPLACEMENT PROBE NO. 20 


Fig. 9 Thrust Bearing Film Thickness Variation With Time Under 
Externally-Imposed Shaped Random Vibrations According 
To NASA Spec 417-2-C-3.5 


MTI-14032 





Fig. 10 Schematic Of BRU Simulator With Static Rotor Force Designation 
And Typical Thrust Runner Position During Vibration Tests 


MTI-14052 



450 



TIME -MINUTES 


Fig. 11 Measured Temperatures In BRU Simulator Components With 
Axial Random Excitation (Shaft Speed 36,000 rpm) 


MTI-12291 
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Fig. 12 Casing-To-Pad Leading Edge Amplitudes For Flex-Mounted Turbine Journal 
Bearing Pad Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 





Fig. 13 Thrust Bearing Film Thickness Variation Under Externally-Imposed 
Shaped Random Vibrations According to NASA Spec 417-2-C-3.5 



INPUT ACC. POWER DENSITY : 0.54 G-RMS INPUT ACC. POWER DENSITY ■ 1.52 G-RMS 



Fig. 14 Amplitude Probability Distribution Of Pad-To-Shaft Pivot Film Thickness 
Variation For Flex-Mounted Compressor Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations According to NASA 
Spec 417-2-C-3.5 



INPU T ACC. POWER DENSITY = 0.54 G-RMS INPUT ACC. POWER DENSITY =1.52 G-RMS 




Fig. 15 Amplitude Probability Distribution Of Pad-To-Shaft Pivot Film Thickness 
Variation For Flex-Mounted Turbine Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations According To 
NASA Spec 417-2-G-3.5 



INPUT ACC. POWER DENSITY ■ 0.54 G-RMS INPUT ACC. POWER DENSITY = 1.52 G-RMS 
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Fig. 16 Amplitude Probability Distribution Of Pad-To-Shaft Pivot Film Thickness 
Variation For Solid-Mounted Turbine Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations According To NASA 

Spec 417-2-C-3.5 MTI-14045 



INPUT ACC. POWER DENSITY = 0.54 6-RMS v INPUT ACC. POWER DENSITY : 1.52 G-RMS 
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Fig. 17 Amplitude Probability Distribution Of Compressor Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 



INPUT ACC. POWER DENSITY 0.54 G-RMS INPUT ACC. POWER DENSITY : 1.52 G-RMS 
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Fig. 18 Amplitude Probability Distribution Of Turbine Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2 -C-3. 5 



INPUT ACC. POWER DENSITY ; 0.54 G-RMS , c INPUT ACC. POWER DENSITY = 1.52 G-RMS 



Fig. 19 Amplitude Probability Distribution Of Casing-To-Pad Leading Edge 
Amplitudes For Flex-Mounted Turbine Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations According To NASA 
Spec 417-2-G-3 .5 



INPUT ACC. POWER DEN SITY • 0.54 6- RMS INPUT ACC. POWER DENSITY ? 1.52 G-RMS 




Fig. 20 Amplitude Probability Distribution Of Compressor Journal Rot< 
Amplitudes (Casing-To-Shaft) Under Externally-Imposed Shaped 
Random Vibrations According To NASA Spec 417-2-C-3.5 



INPUT ACC. POWER DENSITY = 0.54 G-RMS , c INPUT ACC. POWER DENSITY = 1.52 G-RMS 
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Fig. 21 Amplitude Probability Distribution Of Turbine Journal Rotor 
Amplitudes (Casing-To-Shaft) Under Externally-Imposed 
Shaped Random Vibrations According To NASA Spec 
417-2-C-3 .5 



INPUT ACC. POWER DENSITY j 0.54 G-RMS , c INPUT ACC. POWER DENSITY = 1.52 G-RMS 





Fig. 22 Amplitude Probability Distribution Of Thrust Bearing Film Thickness 
Variation Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 



INPUT ACC. POWER DENSITY = 0.54 6-RMS , „ INPUT ACC. POWER DENSITY = 1.52 6- RMS 




Fig. 23 Amplitude Probability Distribution Of Thrust Bearing Gimbal 
Amplitudes Under Externally-Imposed Shaped Random 
Vibrations According To NASA Spec 417-2-C-3.5 



INPUT ACC. POWER DENSITY = 0.54 G-RMS INPUT ACC. POWER DENSITY : 1.52 G-RMS 
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INPUT ACC. POWER DENSITY = 0.54 G-RMS ir INPUT ACC. POWER DENSITY : 1.52 6- RMS 




Fig. 25 Amplitude Probability Distribution Of Pad-To-Shaft Pivot Film 
Thickness Variation For Flex-Mounted Turbine Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 
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Fig. 26 Amplitude Probability Distribution Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Compressor Journal 
Bearing Pad Under Externally-Imposed Shaped Random 
Vibrations According To NASA Spec 417-2-C-3.5 
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Fig. 27 Amplitude Probability Distribution Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Turbine Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 
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Fig. 28 Amplitude Probability Distribution Of Compressor Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 
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Fig. 29 Amplitude Probability Distribution Of Turbine Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 
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Fig. 30 Amplitude Probability Distribution Of Casing-To-Pad Leading Edge 
Amplitudes For Flex-Mounted Turbine Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations According To NASA 
Spec 417-2-C-3 .5 -• ' 
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Fig. 31 Amplitude Probability Distribution Of Compressor Journal Rotor 
Amplitudes (Casing-To-Shaft) Under Externally-Imposed Shaped 
Random Vibrations According To NASA Spec 417-2-C-3.5 





Fig. 32 Amplitude Probability Distribution Of Turbine Journal Rotor 
Amplitudes (Casing-To-Shaf t) ' Under Shaped Random Vibrations 
According To NASA Spec 417 -2-C-3.5 
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Fig. 33 Amplitude Probability Distribution Of Thrust Bearing Film Thickness 
Variation Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 
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Fig. 34 Amplitude Probability Distribution Of Thrust Bearing Gimbal 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
According To NASA Spec 417-2-C-3.5 
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Fig. 36 Random Vibration Displacement Power Spectral Density Test Specification 
417-2 (Rev. C) For Electrical Generating System Components (Operating) 
During Space Flight Operation 
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Fig. 37 Frequency Distributions Obtained With 6, 1.5 And 0.3 Hz Bandwidth For 64 And 
Ensembles Of Pad-To-Shaft Film Thickness Variation Measurements For The Solid- 
Mounted Compressor Journal Bearing Pad Under Externally-Imposed Shaped Random 
Vibrations (1.52 g rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 39 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Flex-Mounted Turbine Journal Bearing 
Pad Under Externally- Imposed Shaped Random Vibrations (1.52 g 
rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 40 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Compressor Journal 
Bearing Pad Under Externally-Imposed Shaped Random Vibra- 
tions (0.54 g rms Input) According To NASA Spec 417-2-C-3.5 




Fig. 41 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Compressor Journal 
Bearing Pad Under Externally- Imposed Shaped Random Vibra- 
tions (1.52 g rms Input) According To NASA Spec 417-2-C-3.5 






ZH/ Z NI-X1ISN3Q IVdlDBdS 83M0d lN3H30V"ldSIQ 


Vd 


£f 

'o 


ro 


ZH / W - A1ISN3Q 1 VH 103 dS H3M0d lN3N30V“ldSia 

V ^ 


Fig. 42 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Turbine Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations (0.54 g 
rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 43 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Turbine Journal Bearing 
\ Pad Under Externally-Imposed Shaped Random Vibrations (1.52 g 

rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 44 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
(0.54 g rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 45 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
(1.52 g rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 46 Frequency Distribution (5 - 2000 Hz) Of Casing-To-Pad Leading Edg 
Amplitudes For Flex-Mounted Turbine Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations (0.54 g rms Input) 
According To NASA Spec 417-2-C-3.5 
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Fig. 47 Frequency Distribution (5 - 2000 Hz) Of Casing-To-Pad Leading Edge 
Amplitudes For Flex-Mounted Turbine Journal Bearing Pad Under 

Externally-Imposed Shaped Random Vibrations (0.75 g rms Input) mti-i3683 

According To NASA Spec 417-2-C-3.5 
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INPUT— VIBRATION 

DISPLACEMENT P.S.D. 
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Fig. 49 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Rotor 
Amplitudes (Casing-To-Shaft ) Under Externally-Imposed Shaped 
Random Vibrations (0.54 g rms Input) According To NASA Spec 
417-2-C-3 .5 
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Fig. 50 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Rotor 
Amplitudes (Casing-To-Shaf t) Under Externally-Imposed Shaped 
Random Vibrations (1.52 g rms Input) According To NASA Spec 
417-2-C-3 .5 
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Fig. 52 Frequency Distribution (5 - 2000 Hz) Of Turbine Journal Rotor 
Amplitudes (Casing-To-Shaf t) Under Externally-Imposed Shaped 
Random Vibrations (1.52 g rms Input) According To NASA Spec 
417-2-C-3 .5 
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Fig. 53 Frequency Distribution (5 - 500 Hz) Of Turbine Journal Rotor 
Amplitudes (Casing-To-Shaf t) Under Externally-Imposed Shaped 
Random Vibrations (1.52 g rms) According To NASA Spec 417-2- 
C-3.5 





Fip. 54 Frequency Distribution (5 - 2000 Hz) Of Thrust Bearing Filir> Thickness 
Variation Under Externally-Imposed Shaped Random Vibrations (0.54 p 
rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 56 Frequency Distribution (5 - 2000 Hz) Of Thrust Bearing Gimbal 
Amplitudes Under Externally- Imposed Shaped Random Vibrations 
(0.54 g rms Input) According To NASA Spec 417-2-C-3.5 
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Fiq. 57 Frequency Distribution (5 - 2000 Hz) Of Thrust Rearing Gimbal 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
(1,52 o rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 58 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film Thicknes 
Variation For Flex-Mounted Compressor Journal Bearing Pad Under 
Externally- Imposed Sinusoidal Vibrations (1.52 g rms Input) 

According To NASA Spec 417-2-C-3.5 
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Fig. 59 Frequency Distribution (5 - 500 Hz) Of Pad-To-Shaft Pivot Film Thickness 
Variation For Flex-Mounted Compressor Journal Bearing Pad Under 
Externally-Imposed Sinusoidal Vibrations (1.52 g rms Input) 

According To NASA Spec 417-2-C-3.5 





ZH/ 2 W -A1ISN3Q im33dS H3M0d iN3W30VldSIQ 


Fig. 60 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Flex-Mounted Turbine Journal Bearing Pad 
Under Externally-Imposed Shaped Random Vibrations (1.35 g rms 
Input) According To NASA Spec 417-2-C-3.5 
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Fig. 61 Frequency Distribution (5 - 500 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Flex-Mounted Turbine Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations (1.35 
g rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 62 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Compressor Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations (0.54 g 
rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 63 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Compressor Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations (1.52 g 
rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 64 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Turbine Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations (0.54 g 
rms Input) According To NASA Spec 417-2-C-3.5 












Fig. 65 Frequency Distribution (5 - 2000 Hz) Of Pad-To-Shaft Pivot Film 
Thickness Variation For Solid-Mounted Turbine Journal Bearing 
Pad Under Externally-Imposed Shaped Random Vibrations (1.35 g 
rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 66 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Flexure 
Amplitudes Under Externally-Imposed Shaped Random Vibrations (0.54 
g rms Input) According To NASA Spec 417-2-C-3.5 
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Fig. 68 Frequency Distribution (5 - 2000 Hz) Of Casing-To-Pad Leading Edg 
Amplitudes For Flex-Mounted Turbine Journal Bearing Pad Under 
Externally-Imposed Shaped Random Vibrations (0.54 g rms Input) 
According To NASA Spec 417-2-C-3.5 
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Fig. 70 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Rotor 
Amplitudes (Casing-To-Shaft) Under Externally-Imposed Shaped 
Random Vibrations (0.54 g rms Input) According To NASA Spec 
417-2-C-3 .5 
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Fig. 71 Frequency Distribution (5 - 2000 Hz) Of Compressor Journal Rotor 
Amplitudes (Casing-To-Shaf t) Under Externally-Imposed Shaped 
Random Vibrations (1.52 g rms Input) According To NASA Spec 

417-2-C-3.5 MTI-13672 





Fig. 72 Frequency Distribution (5 - 500 Hz) Of Compressor Journal Rotor 
Amplitudes (Casing-To-Shaft) Under Externally-Imposed Shaped 
Random Vibrations (1.52 g rms Input) According To NASA Spec 
417-2-C-3 .5 





Fig. 73 Frequency Distribution (5 - 2000 Hz) Of Thrust Bearing Film Thickness 
Variation Under Externally-Imposed Shaped Random Vibrations (0.54 g 
rms Input) According To NASA Spec 417-2-C-3.5 
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uency Distribution (5 - 2000 Hz) Of Thrust Bearing Film Thickness 
ation Under Externally-Imposed Shaped Random Vibrations (1.35 g 
Input) According To NASA Spec 417-2-C-3.5 
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Fig. 75 Frequency Distribution (5 - 2000 Hz) Of Thrust Bearing Gimbal 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
(1.35 g rms Input) According To NASA Spec 417-2-C-3.5 



* L 

'© 


'2 


*2 


% 


TSr 

'o 


ZH fJH - UISN3Q 1VM133d$ M3M0d iN3N33V3dSIQ 


Fig. 76 Frequency Distribution (5 - 500 Hz) Of Thrust Bearing Gimbal 
Amplitudes Under Externally-Imposed Shaped Random Vibrations 
(1.35 g rms Input) According To NASA Spec 417-2-C-3.5 





Fig. 77 Lumped-Mass Model Used For Calculation Of Axial Vibration 
Response Of BRU Rotor-Bearing System To Specified Random 
Base Excitation 
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Fig. 79 Random Vibration Power Spectral Density Test Specification 417-2 
(Rev. C) For Electrical Generating System Components (Operating) 
And Actual BRU Test Levels 


MTI-14053 
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